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ABSTRACT 

We investigate the impact of both slow and fast polarization modulation strategies on the 
science return of upcoming ground-based experiments aimed at measuring the _B-mode po- 
larization of the cosmic microwave background. Using detailed simulations of the qover 
experiment, we compare the ability of modulated and un-modulated observations to recover 
the signature of gravitational waves in the polarized CMB sky in the presence of a number 
of anticipated systematic effects. The general expectations that fast modulation is helpful in 
mitigating low-frequency detector noise, and that the additional redundancy in the projection 
of the instrument's polarization sensitivity directions onto the sky when modulating reduces 
the impact of instrumental polarization, particularly for fast modulation, are borne out by our 
simulations. Neither low-frequency polarized atmospheric fluctuations nor systematic errors 
in the polarization sensitivity directions are mitigated by modulation. Additionally, we find no 
significant reduction in the effect of pointing errors by modulation. For a QOVER-like experi- 
ment, pointing jitter should be negligible but any systematic mis-calibration of the polarization 
coordinate reference system results in significant E-B mixing on all angular scales and will 
require careful control. We also stress the importance of combining data from multiple detec- 
tors in order to remove the effects of common-mode systematics (such as un-polarized 1// 
atmospheric noise) on the measured polarization signal. Finally we compare the performance 
of our simulated experiment with the predicted performance from a Fisher analysis. We find 
good agreement between the (optimal) Fisher predictions and the simulated experiment ex- 
cept for the very largest scales where the power spectrum estimator we have used introduces 
additional variance to the S-mode signal recovered from our simulations. In terms of detect- 
ing the total B-mode signal, including lensing, the Fisher analysis and the simulations are in 
excellent agreement. For a detection of the primordial B-mode signal only, using an input 
tensor-to-scalar ratio of r = 0.026, the Fisher analysis predictions are ~ 20 per cent, better 
than the simulated performance. 

Key words: methods: statistical - methods: analytical - cosmology: theory - cosmic mi- 
crowave background 



1 INTRODUCTION 

There is currently a great deal of interest in the rapidly evolving 
field of observations of the polarization of the cosmic microwave 
background (CMB). This interest stems from the fact that such ob- 
servations have the potential to discriminate between inflation and 
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other early-universe models through their ability to constrain an 
odd-parity 5-mode polarization component induced by a stochas- 
tic background of gravitational waves at the time of last scatter- 
ing (Kamionkowsky, Kosowksy & Stebbins 1997; Seljak & Zal- 
darriaga 1997). 

From an observational point of view, we are still a long way 
off from a detection of this B-mode signature of inflation. However, 
much progress has been made recently with the detection of the 
much stronger _E-mode polarization signal on large scales by the 
wmap experiment (Page et al. 2007; Nolta et al. 2008). On smaller 
scales, a growing number of balloon-borne and ground-based ex- 
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periments have also measured B-mode polarization including dasi 
(Leitch et al. 2005), CBI (Sievers et al. 2007), boomerang (Mon- 
troy et al. 2006), maxipol (Wu et al. 2007), CAPMAP (Bischoff et 
al. 2008) and QUaD (Pryke et al. 2009). Most recently, the high 
precision measurement of small-scale polarization by the QUaD ex- 
periment has, for the first time, revealed a characteristic series of 
acoustic peaks in the B-mode spectrum and put the strongest up- 
per limits to date on the small-scale B-mode polarization signal 
expected from gravitational lensing by large-scale structure. 

Building on the experience gained from these pioneering ex- 
periments, a new generation of experiments is now under construc- 
tion with the ambitious goal of observing the primordial B-mode 
signal. Observing this signal is one of the most challenging goals 
of modern observational cosmology. There are a number of rea- 
sons why these types of observations are so difficult. First and fore- 
most, the sought-after signal is expected to be extremely small - 
in terms of the tensor-to-scalar ratio 1 , the RMS polarization sig- 
nal from primordial B-modes is 0.4-/F /xK, and the current 95 per 
cent limit r < 0.22 from wmap temperature and B-mode po- 
larization plus distance indicators (Komatsu et al. 2009) implies 
an RMS < 180 nK. Secondly, polarized emission from our own 
galaxy and from extra-Galactic objects act as a foreground contam- 
inant in observing the CMB polarized sky. Although our knowledge 
of such polarized foregrounds is currently limited, particularly at 
the higher frequencies > 100 GHz of relevance to bolometer ex- 
periments, models suggest that such contamination could be an or- 
der of magnitude larger than the sought-after signal on the largest 
scales (e.g. Amblard, Cooray & Kaplinghat 2007). Thirdly, gravi- 
tational lensing by large-scale structure converts B-modes into B- 
modes on small to medium scales (see Lewis & Challinor 2006 for 
a review) and acts as a source of confusion in attempts to mea- 
sure the primordial B-mode signal (Knox & Song 2002; Kesden, 
Cooray & Kamionkowski 2002). Note however that the lensing B- 
mode signal is a valuable source of cosmological information in its 
own right and can be used to put unique constraints on dark en- 
ergy and massive neutrinos (e.g. Kaplinghat, Knox & Song 2003; 
Smith, Challinor & Rocha 2006). Unfortunately these latter two 
effects (foreground contamination and weak gravitational lensing) 
contrive in such a way as to render B-mode polarization observa- 
tions subject to contamination on all angular scales (the primordial 
B-mode signal is dominated by foregrounds on large scales whilst 
on smaller scales it is swamped by the lensing signal). Last, but 
not least, exquisite control of systematic and instrumental effects 
will be required, to much better than 100 nK, before any detection 
of B-modes can be claimed. The sought-after signal is so small 
that systematic and instrumental effects considered negligible for 
an exquisitely precise measurement of B-modes say, could poten- 
tially ruin a detection of B-modes, if left uncorrected. One possible 
approach to mitigating some of these systematics in hardware is to 
modulate the incoming polarization signal such that it is shifted to 
higher frequency and thus away from low-frequency systematics 
which would otherwise contaminate it. There are a number of tech- 
niques for achieving this including the use of a rotating half-wave 



1 Our normalisation conventions follow those adopted in the CAMB 
code (Lewis, Challinor & Lasenby 2000), so that r is the ratio of primordial 
power spectra for gravitational waves and curvature perturbations. Explic- 
itly, for slow-roll inflation in a potential V(<j>), 2r fs M^JV" (4>)/V((f>)] 2 
where M Pl = 2.436 X 10 18 GeV/c 2 is the reduced Planck mass. 



plate (HWP; e.g. Johnson et al. 2007), phase-switching, or Faraday 
rotation modulators (Keating et al. 2003). 

In this paper, we investigate the ability of modulation tech- 
niques that are either slow or fast with respect to the temporal vari- 
ation of the signal to mitigate a range of possible systematic effects. 
Our analysis is based on simulations, and the subsequent analysis, 
of data from a ground-based CMB B-mode polarization experi- 
ment. Previous investigations of the impact of systematic effects on 
B-mode observations include the analytic works of Hu, Hedman & 
Zaldariagga (2003), O'Dea, Challinor & Johnson (2007) and Shi- 
mon et al. (2008), as well as the simulation-based analysis of Mac- 
Tavish et al. (2008) who based their study on signal-only simula- 
tions of the spider experiment. The simulation work presented here 
is complementary to these previous analyses but we also take our 
analysis further by including realistic noise in our simulations — 
we are thus able to quantify not only any bias found, but also any 
degradation of performance due to the presence of systematic ef- 
fects. Our work makes use of a detailed simulation pipeline which 
we have created in the context of the QOVER experiment. Although 
the precise details of our simulations are specific to QOVER, our 
general conclusions regarding the impact of modulation on a vari- 
ety of systematic effects are relevant to all upcoming ground-based 
B-mode experiments and many of them are also relevant for both 
balloon-borne and space-based missions. 

The paper is organised as follows. In Section 2, we review 
the relevant upcoming B-mode polarization experiments. Section 
3 describes our simulation technique and the systematic effects 
we have considered. In Section 4, we describe the map-making 
and power spectrum estimation techniques that we use to analyse 
the simulated data. Section 5 presents the results from our main 
analysis of systematic effects. We discuss our results in Section 6 
where, for clarity, we group the possible systematic effects con- 
sidered into those that are, and those that are not, mitigated by a 
modulation scheme. In this section, we also demonstrate the im- 
portance of combining information from multiple detectors during 
analysis and compare the simulated performance of Clover to the 
predicted performance from a Fisher analysis. Our conclusions are 
summarised in Section 7. Finally, Appendix A develops a simple 
model of the spurious B-mode power produced by pointing jitter 
for experiments with a highly redundant scan strategy. 



2 CMB POLARIZATION EXPERIMENTS 

A host of experiments are currently under construction (one is, in 
fact, already observing) with their primary goal being to constrain 
the tensor-to-scalar ratio, and thus the energy scale of inflation, 
through observations of the B-mode component of the CMB. Here, 
we give a short summary of the planned experiments. 

( i) Ground-based experiments: 

• bicep/bicep-2/keck array: The bicep experiment (Yoon et al. 
2006) has recently completed its third and final season of opera- 
tion from the South Pole. The experiment consisted of a total of 98 
polarization-sensitive bolometers (PSBs) at 100 and 150 GHz. The 
optical design is very clean but with the downside of poor resolu- 
tion - 45 arcmin full- width at half-maximum (FWHM) at 150 GHz 
- which limits the target multipole range to t < 300. The stated 
target sensitivity is r = 0.1. In the first observing season, three 
100 GHz pixels and three 150 GHz pixels were equipped with 
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Faraday rotation modulators and two pixels operated at 220 GHz. 
bicep-2 will consist of an upgrade to the bicep telescope with a 
150 GHz 512-element array of antenna-coupled detectors. It will 
be deployed to the South Pole in November, 2009. The KECK ar- 
ray will consist of three BlCEP-2-like telescopes (at 100, 150 and 
220 GHz). It is hoped to be installed on the dasi mount (previ- 
ously occupied by QUaD) in November 2010. The nominal goal of 
this array is r = 0.01. 

• QOVER: For an up-to-date overview see North et al. (2008). 
QOVER is a three-frequency (97, 150 and 225 GHz) instrument 
to be sited at Pampa La Bola in the Atacama desert in Chile. It 
will have 576 single-polarization transition-edge sensors (TES), 
split equally among the three frequencies. The beam size of ~ 
5.5 arcmin FWHM at 150 GHz will sample the multipole range 
25 < I < 2000. The target sensitivity is r ~ 0.03 and the polariza- 
tion signal will be modulated with a HWP The 97 GHz instrument 
is expected to be deployed to Chile in late 2009 with the combined 
150/225 GHz instrument to follow soon after in 2010. 

• quiet: See Samtleben (2008) for a recent overview, quiet is 
unique among planned S-mode experiments in that it uses pseudo- 
correlation HEMT-based receivers rather than bolometers. It will 
observe from Chile using the CBI mount - a planned second phase 
will involve upgrading to ~ 1000 element arrays and relocation of 
the 7-m Crawford Hill antenna from New Jersey to Chile, quiet 
will observe at 40 and 90 GHz. The target sensitivity for the second 
phase is r ~ 0.01. 

• POLARBEAR: A three-frequency (90, 150 and 220 GHz) single- 
dish instrument to be sited in the Inyo Mountains, CA for its first 
year of operation, after which it will be relocated to the Atacama 
desert in Chile. It will use 1280 TES bolometers at each frequency 
with polarization modulation from a HWP. The planned beam size 
is 4 arcmin (FWHM) at 150 GHz. The target sensitivity is r = 
0.015 for the full instrument. 

• brain: See Charlassier (2008) for a recent review, brain is a 
unique bolometric interferometer project (c.f. dasi, CBl) to be sited 
on the Dome-C site in Antarctica. The final instrument will have 
~ 1000 bolometers observing at 90, 150 and 220 GHz. brain will 
be primarily sensitive to multipoles 50 < t < 200. The full ex- 
periment is planned to be operational in 201 1 and the stated target 
sensitivity is r = 0.01. 

( ii) Balloon-borne experiments 

• ebex: See Oxley et al. (2004) for a summary, ebex will ob- 
serve at 150, 250 and 410 GHz and will fly a total of 1406 TES with 
HWP modulation. The angular resolution is 8 arcmin and the target 
multipole range is 20 < I < 2000. The stated target sensitivity 
is r — 0.02. A test flight is planned for 2009 and a long-duration 
balloon (LDB) flight is expected soon after. 

• spider: See Crill et al. (2008) for a recent description, spider 
will deploy ~ 3000 antenna-coupled TES observing at 96, 145, 
225 and 275 GHz, with a beam size of ~ 40 arcmin at 145 GHz. 
The target multipole range is 10 < I < 300. A 2-6 day first flight 
is planned for 2010. The target sensitivity is r — 0.01. Signal mod- 
ulation will be provided by a (slow) stepped HWP and fast gondola 
rotation. 



• piper: This balloon experiment will deploy a focal plane of 
5120 TES bolometers in a backshort-under-grid (BUG) configu- 
ration. Each flight of piper will observe at a different frequency, 
covering 200, 270, 350 and 600 GHz after the four planned flights. 
The beam size is ~ 15 arcmin, corresponding to a target multipole 
range £ < 800. The first element of the optical system is a variable 
polarization modulator (VPM). The entire optical chain, including 
the modulators, are cooled to 1.5 K so that piper observes at the 
background limit for balloon altitudes. Including removal of fore- 
grounds, the experiment has the sensitivity to make a 2a detection 
of r = 0.007. The first flight is scheduled for 2013. 

( Hi) Space missions 

• planck: See the publication of the Planck Collaboration 
(2006) for a detailed review of the science programme, planck 
will measure the temperature in nine frequency bands, seven of 
which will have (some) polarization sensitivity. The polarized 
channels (100, 143, 217 and 353 GHz) of the high-frequency in- 
strument (HFI) use similar PSBs to those deployed on boomerang 
and have beam sizes 5-9.5 arcmin. For low r, sensitivity to pri- 
mordial gravitational waves will mostly come from the large-angle 
reionisation signature (Zaldarriaga 1997) and r = 0.05 may be 
possible if foregrounds allow, planck will be sensitive to the mul- 
tipole range 2 < I < 3000 and is scheduled to launch in 2009. The 
HFI has no active or fast signal modulation (i.e. other than scan- 
ning). 

• CMBPOL/bpol: Design studies have been conducted for satel- 
lite mission(s) dedicated to measuring primordial B-modes com- 
prising ~ 2000 detectors with the ability to measure 0.001 < r < 
0.01 if foregrounds allow it (Bock et al. 2008; de Bernardis et al. 
2008). The timescale for launch of any selected mission is likely 
beyond 2020. 



3 SIMULATIONS 

The QOVER experiment will consist of two telescopes - a low 
frequency instrument with a focal plane consisting of 192 single- 
polarization 97 GHz TES detectors, and a high frequency instru- 
ment with a combined focal plane of 150 and 225 GHz detectors 
(192 of each). Note that we have not included foreground contami- 
nation in our simulations, so for this analysis, we consider only the 
150 GHz detector complement - the corresponding reduction in 
sensitivity will approximate the effect of using the multi-frequency 
observations to remove foregrounds. Figure 1 shows the arrange- 
ment of the 150 and 225 GHz detectors on the high frequency fo- 
cal plane. The detectors are arranged in detector blocks consist- 
ing of eight pixels each. Each pixel consists of two TES detectors 
which are sensitive to orthogonal linear polarizations. The polariza- 
tion sensitivity of the eight detector pairs within a block are along, 
and at right angles to, the major axis of their parent block. The de- 
tector complement, both at 150 and 225 GHz, therefore consists of 
three 'flavours' of pixels with different polarization sensitivity di- 
rections. The 97 GHz focal plane (not shown) has a similar mix of 
detector orientations. 

3.1 Simulation parameters 

Because simulating the full QOVER experiment is computationally 
demanding (a single simulation of a two-year campaign at the full 
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Figure 1. Layout of detectors on the QOVER high frequency focal plane. 
Each point indicates a pixel comprising two TES detectors sensitive to or- 
thogonal linear polarizations. The polarization sensitivity directions of de- 
tectors within each block are along, and at right angles to, the major axis of 
the block. The outlined 150 GHz detector blocks at the centre of the array 
are used in the simulations described here. The field of view of the entire 
array is ~ 5 deg 2 . 



Clover data rate would require ~ 10 4 CPU-hours), we have scaled 
some of the simulation parameters in order to make our analysis 
feasible. 

(i) We simulate only half the 150 GHz detector complement and 
have scaled the noise accordingly. We have verified with a restricted 
number of simulations using all the 150 GHz detectors that the 
marginally more even coverage obtained across our field using all 
detectors has little or no impact on our results or conclusions. The 
150 GHz detector blocks used in our simulations are indicated in 
Fig. 1 and include all three possible orientations of pixels on the 
focal plane. 

(ii) The QOVER detectors will have response times of ~ 200 /is 
and so the data will be sampled at ~ 1 kHz in order to sample 
the detector response function adequately. Simulating at this rate 
is prohibitive so we simulate at a reduced data rate of 100 Hz. For 
the QOVER beam size (FWHM = 5.5 arcmin at 150 GHz) and out- 
chosen scan speed (0.25°/s), this data rate is still fast enough to 
sample the sky signal adequately. 

(iii) QOVER will observe four widely separated fields on the sky, 
each covering an area of ~ 300 deg 2 , over the course of two years. 
Two of the fields are in the southern sky and two lie along the equa- 
tor. For our analysis, we observe each of the four fields for a single 
night only and, again, we have scaled the noise levels to those ap- 
propriate for the full two-year observing campaign. 



3.2 Observing strategy 

Although optimisation of the observing strategy is not the focus of 
this work, a number of possible strategies have been investigated 
by the Clover team. For our analysis, we use the most favoured 



scan strategy at the time of writing. To minimise rapid variations in 
atmospheric noise, the two telescopes will scan back and forth in 
azimuth at constant elevation, allowing the field to rise through the 
chosen observing elevation. Every few hours, the elevation angle 
of the telescopes will be re-pointed to allow for field-tracking. Al- 
though the precise details of the scan are likely to change, the gen- 
eral characteristics of the scan and resulting field coverage prop- 
erties will remain approximately the same due to the limitations 
imposed by constant-elevation scanning and observing from Ata- 
cama. The Clover telescopes are designed with the capability of 
scanning at up to 10°/s. However, for our analysis where we have 
considered QOVER operating with a rotating HWP (Section 3.7), 
we have chosen a relatively slow scan speed of 0.25°/s in light of 
the HWP rotation frequency which we have employed (fx = 3 Hz). 
Although the mode of operation of a HWP on QOVER is still under 
development, a continuously rotating HWP is likely to be restricted 
to rotation frequencies of fx < 5 Hz due to mechanical constraints 
(with current cryogenic rotation technologies, fast rotation, > 5 Hz, 
could possibly result in excessive heat generation). Figure 2 shows 
the coverage maps for a single day's observing on one of the south- 
ern fields and on one of the equatorial fields. The corresponding 
maps for the other two fields are broadly similar. Note that, in the 
real experiment, we expect that somewhat more uniform field cov- 
erage than that shown in Fig. 2 will be achievable by employing 
slightly different scan patterns on different days. 



3.3 Signal simulations 

We generate model TT, EE, TE and BB CMB power spectra 
using CAMB (Lewis, Challinor & Lasenby 2000). The input cos- 
mology used consists of the best-fit standard ACDM model to the 
5-year wmap data set (Hinshaw et al. 2009), but with a tensor- 
to-scalar ratio of r = 0.026, chosen to match the QOVER 'target' 
value. Realisations of CMB skies from these power spectra are then 
created using a modified version of the healpix 2 software (Gorski 
et al. 2005). Our simulations include weak gravitational lensing 
but ignore its non-Gaussian aspects. Using only Gaussian simu- 
lations means that we slightly mis-estimate the covariance matrices 
of our power spectrum estimates, particularly the _B-mode covari- 
ances (Smith, Hu & Kaplinghat 2004; Smith, Challinor & Rocha 
2006). For the QOVER noise levels, this is expected to have a neg- 
ligible impact on the significance level of the total B-mode signal. 

As part of the simulation process, the input CMB signal is con- 
volved with a perfect Gaussian beam with FWHM = 5.5 arcmin. 
Note that an important class of systematic effects which we do not 
consider in this paper are those caused by imperfect optics; see the 
discussion in Section 3.6. 

For our analysis of the simulated data sets (Section 4) we have 
chosen to reconstruct maps of the Stokes parameters with a map 
resolution of 3.4 arcmin (healpix 2V s ide = 1024). Note that this 
pixel size does not fully sample the beam and it is likely that we will 
adopt Af s ido = 2048 for the analysis of real data. In order to iso- 
late the effects of the various systematics we have considered, our 
simulated CMB skies have therefore been created at this same res- 
olution — we can then be sure that any bias found in the recovered 
CMB signals is due to the systematic effect under consideration 



See http://healpix.jpl.nasa.gov 
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Figure 2. Hit-maps for one of the southern fields (left; RA 09:30 hrs , Dec -40.00°) and one of the equatorial fields (right; RA 04:00 hrs, Dec 0.00°) for a 
single day's observation with half the 150 GHz detector complement. The central part of the fields (shown in yellow and red) are roughly 20° in diameter. 
These maps have been constructed using a HEALPIX resolution of N s ^ c = 1024 corresponding to a pixel size ~3.4 arcmin. 



rather than due to a poor choice of map resolution 3 . Note that our 
adopted procedure of simulating and map-making at identical res- 
olutions, although useful for the specific aims of this paper, is not 
a true representation of a CMB observation. For real observations, 
pixelisation of the CMB maps will introduce a bias to the measured 
signal on scales comparable to the pixel size adopted. 

Using the pointing registers as provided by the scan strategy 
and after applying the appropriate focal plane offsets for each de- 
tector, we create simulated time-streams according to 



di = [T(6) + Q(0) cos(2^) + U(0) sin(20i)] A 



(1) 



where 6 denotes the pointing and T, Q and U are the sky signals 
as interpolated from the input CMB sky map. The polarization an- 
gle, fa is, in general, a combination of the polarization sensitivity 
direction of each detector, any rotation of the telescope around its 
boresight, the direction of travel of the telescope in RA-Dec space 
and the orientation of the half-wave plate, if present. 



3.4 Noise simulations 

The QOVER data will be subject to several different noise sources. 
Firstly, photon loading from the telescope, the atmosphere and the 
CMB itself will subject the data to uncorrelated random Gaus- 
sian noise. Secondly, the TES detectors used in QOVER are sub- 
ject to their own sources of noise which will possibly include 
low-frequency 1/f behaviour and correlations between detectors. 
Thirdly, the atmosphere also has a very strong 1/f component 
which will be heavily correlated across the detector array. Fortu- 
nately, the 1/f component of the atmosphere is known to be al- 
most completely un-polarized and so can be removed from the po- 
larization analysis by combining data from multiple detectors (see 
Section 6.2). 



3 We have verified that spurious B-modes generated through the pixelisa- 
tion are negligible for the QOVER noise levels. 



The white-noise levels due to loading from the instrument, at- 
mosphere and CMB have been carefully modelled for the case of 
QOVER observations from Atacama. We will not present the de- 
tails here, but for realistic observing conditions and scanning ele- 
vations, we have calculated the expected noise-equivalent tempera- 
ture (NET) due to photon noise alone to be « 146 fiK^/s. We add 
this white noise component to our simulated signal time-streams 
for each detector as 

NET 

di — s- di H \J /sampgi, (2) 

where / samp is the sampling frequency and gi is a Gaussian random 
number with fi = and a = 1. Note that the white-noise level in 
the detector time streams is NET/2 since the QOVER detectors are 
half-power detectors (equation 1). 

Using instrument parameters appropriate for the Clover de- 
tectors, we use the small-signal TES model of Irwin & Hilton 
(2005) to create a model noise power spectrum for the detector 
noise. This model includes both a contribution from the super- 
conducting SQUIDs which will be used to record the detector sig- 
nals (e.g. Reintsema et al. 2003) and a contribution from aliasing in 
the Multi Channel Electronics (MCE; Battistelli et al. 2008) which 
will be used to read out the signals. For the instrument parameters 
we have chosen, the effective NET of the detector noise in our sim- 
ulations is approximately equal to the total combined photon noise 
contribution from the atmosphere, the instrument and the CMB. 
Note however that for the final instrument, it is hoped that the de- 
tector NET can be reduced to half that of the total photon noise, 
thus making QOVER limited by irreducible photon loading. The 
Irwin & Hilton (2005) small-signal TES model does not include 
a 1// component to the detector noise so in order to investigate 
the impact of modulation on possible low-frequency detector noise, 
we add a heuristically chosen 1 // component to the detector noise 
model with knee frequencies in the range, 0.01 < /knee < 0.1 Hz. 
The MCE system which will be used to read out the QOVER data 
should have low cross-talk between different channels. However, 
correlations will be present at some level and so we include 10 



© 2009 RAS, MNRAS 000, 1-24 



6 M.L. Brown et al. 



per cent correlations between all of our simulated detector noise 
time-streams. Generally, to simulate stationary noise that is corre- 
lated in time and across the iVdct detectors we proceed as follows. 
Let the noise cross-power spectrum between detector d and d' be 
Pd,d> (/) • Taking the Cholesky decomposition of this matrix at each 
frequency, L dd i (/), defined by 



p d ,d'(/) = X)W'(/)£<i',d''(/)> 



(3) 



we apply L to iVdet independent, white-noise time streams g<i{f) 
in Fourier space, 



9d(f) -^^2L d:d ,(f)g d ,(f). 



(4) 



The resulting time-streams, transformed to real space then pos- 
sess the desired correlations between detectors. Here, we assume 
that the correlations are independent of frequency so that the noise 
cross-power spectrum takes the form 



Pd,d>(f) = Cd,d>P(f), 



(5) 



where the correlation matrix C d>d i = 1 for d = d' and C d>d i =0.1 
otherwise. In practice, we use discrete Fourier transforms to synthe- 
sise noise with periodic boundary conditions (and hence circulant 
time-time correlations). 

We use the same technique to simulate the correlated 1/f 
component of the atmosphere. We have measured the noise proper- 
ties of the atmosphere from data from the QUaD experiment (Hin- 
derks et al. 2009). The 150 GHz frequency channel of QUaD is obvi- 
ously well matched to the QOVER 150 GHz channel although QUaD 
observed the CMB from the South Pole rather than from Atacama. 
Although there are significant differences between the properties of 
the atmosphere at the South Pole and at Atacama (e.g. Bussmann, 
Holzapfel & Kuo 2005), the QUaD observations still represent the 
best estimate of the 1/f noise properties of the atmosphere avail- 
able at present. A rough fit of the QUaD data to the model, 



P(f) = NET 2 



1 + 



A 



/ 



(6) 



yields a knee frequency, /knee = 0.45 Hz and spectral index, 
a — 2.5. Using this model power spectrum we simulated 1/ / at- 
mospheric noise correlated across the array in exactly the same way 
as was used for the detector noise. Fortunately, the 1/f compo- 
nent in the atmosphere is almost completely un-polarized. If there 
were no instrumental polarization, detectors within the same pixel 
(which always look in the same direction) would therefore be com- 
pletely correlated with one another and detectors from different 
pixels would also be heavily correlated. For the correlated atmo- 
sphere, we therefore use a correlation matrix given by C d:d / = 1.0 
for \d— d'\ ^ 1 (i.e. for detector pairs) and C d>d ' = 0.5 otherwise. 
In the following sections, as one of the systematics we have investi- 
gated, we relax the assumption that the atmosphere is un-polarized. 

Figure 3 compares the photon, atmospheric 1/f and detector 
noise contributions to our simulated data in frequency space. At low 
frequencies, the noise is completely dominated by the atmospheric 
1/f while the white-noise contributions from photon loading (in- 
cluding the uncorrelated component of the atmosphere) and detec- 
tor noise are approximately equal. Note that for observations with- 
out active modulation, and for the scan speed and observing eleva- 
tions which we have adopted, the "science band" for the multipole 
range 20 < I < 2000 corresponds roughly to 0.01 < / < 1Hz 
in time-stream frequency. In contrast, for our simulations which in- 




Detector noise 
1/f atmospheric noise 
Total photon noise (incl. atms.) 
Total noise 



k (Hz) 

Figure 3. Frequency space comparison between the different noise sources 
in the simulations. The grey line shows the detector noise power spectrum 
(here with a 1// component with knee frequency, /k nC o =0.1 Hz). The 
correlated component to the atmosphere is shown as the dashed line and the 
total photon noise (including atmospheric loading) is shown as the dotted 
line. For the simulation parameters we have adopted, the temperature sky- 
signal from multipoles 20 < i < 2000 appears in the time-stream in the 
frequency range 0.01 < / < 1 Hz. In the absence of fast modulation, 
the polarized sky signal also appears in this frequency range whereas in our 
simulations including fast modulation, the polarized sky signal appears in a 
narrow band centred on 12 Hz. 



elude a continuously rotating HWP, the temperature signal remains 
within the 0.01 < / < 1Hz frequency range but the polarized sky 
signal is moved to a narrow band centred on ~12 Hz, well away 
from both the detector and atmospheric 1/f noise components (see 
Section 3.7 and Fig. 5). Note also that although the atmospheric 
1/f dominates the detector 1// at low frequency, the atmosphere 
is heavily correlated across detectors and can therefore be removed 
by combining detectors (e.g. differencing detectors within a pixel) 
but this is not true for the detector noise which is only weakly cor- 
related between detectors. We demonstrate this in Fig. 4 where we 
plot a five-minute sample of simulated atmospheric and detector 
noise for the two constituent detectors within a pixel. Including the 
atmospheric 1/f component, the effective total NET per detector 
measured from our simulated time-streams is 293 /iK^/sec whilst 
excluding atmospheric 1/f, we measure 210/iK^/s. 



3.5 Detector response 

The Irwin & Hilton (2005) TES small-signal model mentioned 
above also provides us with an estimate of the detector response 
(the conversion from incident power to resultant current in the de- 
tectors). In this model, the power-to-current responsivity, s/(o>), is 
given by 



si(ui) oc 



1 - t+/ti 1 — t-I Tj 

1 + iiOT + 1 + IUJT- 



(7) 



where u) = 27r/ is the angular frequency and r+ and r_ are the 
"rise time" and "fall time" (relaxation to steady state) after a delta- 
function temperature impulse. Here, n is the current-biased ther- 
mal time constant. The impulse-response in the time domain is 



Sl(t) oc 



-*/t+ 



-t/T 



-e(t), 



(8) 
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Figure 4. A five-minute sample of simulated noise time-stream for the two 
detectors within a pixel (denoted "A" and "B") for both the atmospheric 
noise simulations (top) and for the detector noise simulations (bottom). The 
1/f component in the atmospheric noise time-streams is 100 per cent corre- 
lated between the A and B detectors and can be removed entirely by differ- 
encing. In contrast the 1/f component in the detector noise (which is much 
weaker and not noticeable on this plot) is only weakly correlated between A 
and B and is not removed by differencing. For comparison, the signal-only 
time-streams are also plotted in the lower panels as the red curves. 



where O(t) is the Heaviside step function. Note that the constant 
of proportionality in equation (7) (which is also predicted by the 
model) is essentially the calibration (gain) of the detectors. The 
simulated time-streams of equation (1) are converted to detector 
time-streams through convolution with this response function. Note 
that the photon noise and correlated atmospheric noise are added to 
the time-stream before convolution (and so are also convolved with 
the response function) while the detector noise is added directly to 
the convolved time-streams. The Clover detectors are designed to 
be extremely fast with time-constants t± < 1 ms. For our simu- 
lations, we have used time-constants predicted by the small-signal 
TES model for QOVER instrument parameters of r+ = 300 /is and 
r- = 322 [is. 

Note that for our chosen scan speed of 0.25°/s, the effect of 
the response function on the signal component in the simulations 
is small in the absence of fast modulation - the QOVER detectors 
are so fast that signal attenuation and phase differences introduced 
by convolution with the detector response only become important 
at high frequencies, beyond the frequency range of the sky signals. 
For our simulations without fast modulation, sky signals from mul- 
tipoles £ ~ 2000 will appear in the time-stream at ~ 1 Hz where the 
amplitude of the normalised response function is effectively unity 
and the associated phase change is —0.2°. For our simulations in- 
cluding fast modulation, it becomes more important to correct the 
time-stream data for the detector response. In this case, the polar- 
ized sky signals (from all multipoles) appear within a narrow band 
centred on 12 Hz in the time-stream. Here, the amplitude of the re- 
sponse function is still very close to unity but the phase change has 
grown to —2.7° . We include a deconvolution step in the analysis of 
all of our simulated data to correct for this effect. Finally, we note 
that for the ~ 10 per cent errors in the time-constants which we 
have considered (see Section 3.6), the resulting mis-estimation of 
the polarization signal will again be small, even for the case of fast 
modulation. 



In principle we should also include the effect of sample in- 
tegration: each discrete observation is an integral of a continuous 
signal over the sample period. Of course, in the case where down- 
sampled data is simulated the sample integration should include 
the effect of the sample averaging. For integration over a down- 
sampled period A, there is an additional phase-preserving filtering 
by sinc(o;A/2) where u> = 2nf. For our scan parameters, the filter 
is negligible (i.e. unity) for unmodulated data while for modulated 
data the filter can be approximated at the frequency Auj\ , where lu\ 
is the angular rotation frequency of the HWP. This acts like a small 
decrease in the polarization efficiency, with the discretised signal at 
the detector being 

di w i |r(<9) +sinc(2Lj A A) 

x [Q{6) cos(20i) + U{0) sin(20O] }• (9) 

For u\ = 2-7T x 3 Hz and A = 10 ms, the effective polarization 
efficiency is 0.98 and this has the effect of raising the noise level in 
the polarization maps by 2 per cent. We do not include this small 
effect in our simulations, but could easily do so. 

3.6 Systematic effects 

In the analysis that follows, we will investigate the impact of sev- 
eral systematic effects on the ability of a QOVER-like experiment 
to recover an input _B-mode signal. For a reference, we use a suite 
of simulations which contain no systematics. This ideal simulation 
contains the input signal, photon noise, 1 // atmospheric noise cor- 
related across the array (but un-polarized) and TES detector noise 
with no additional correlated 1// component. Additionally, for our 
reference simulation all pointing registers and detector polariza- 
tion sensitivity angles use the nominal values and the signal is con- 
volved with the detector response function using the nominal time 
constants. 

We then perform additional sets of simulations with the fol- 
lowing systematic effects included in isolation: 

• 1// detector noise. We have considered an additional corre- 
lated 1/f component to the detector noise with 1/f knee frequen- 
cies of 0.1, 0.05 and 0.01 Hz. 

• Polarized atmosphere. In addition to the un-polarized atmo- 
sphere present in the reference simulation, we consider a polar- 
ized 1/f component in the atmosphere. To simulate polarized at- 
mosphere, we proceed as described in Section 3.4 but now we add 
correlated 1/f atmospheric noise to the Q and U sky signal time- 
streams such that equation (1) becomes 

d t = X - [T+ (Q + Qf ms ) cos(2^) 

+ (U + Uf ms ) sin(2&)] • (10) 

We take the Q atms and f/ atms atmospheric signals to have the same 
power spectrum as the common-mode atmosphere (equation 6) but 
a factor ten smaller in magnitude. 

• Detector gain errors. We consider three types of gain errors: 
(i) random errors in the gain that are constant in time, uncorrelated 
between detectors and have a 1 per cent RMS; (ii) gain drifts in 
each detector corresponding to a 1 per cent drift over the course of 
a two-hour observation - the start and end gains for each detector 
are randomly distributed about the nominal gain value with an RMS 
of 1 per cent; (iii) systematic A/B gain mis-matches (1 per cent mis- 
match) between the two detectors within each pixel. For this latter 
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systematic, we have applied a constant 1 per cent A/B mis-match to 
all pixels on the focal plane but the direction of the mis-match (that 
is, whether the gain of A is greater or smaller than B) is chosen 
randomly. 

• Mis-estimated polarization sensitivity directions. Random er- 
rors uncorrelated between detectors (including those with the same 
feedhorn) with RMS 0.5° and which are constant in time, and a 
systematic mis-estimation of the instrument polarization coordinate 
reference system by 0.5° are considered. 

• Mis-estimated half-wave plate angles. For the case where we 
consider an experiment which includes polarization modulation 
with a half-wave plate (see Section 3.7), we also consider random 
errors (with RMS 0.5°) in the recorded HWP angle which we ap- 
ply to each 100 Hz time-sample. In addition, we consider a 0.5° 
systematic offset in the half-wave plate angle measurements. 

• Mis-estimated time-constants. The analysis that follows in- 
cludes a deconvolution step to undo the response function of the 
detectors and return the deconvolved sky signal. In all cases, we 
use the nominal time-constant values of r+ = 300 /is and t_ = 
322 /is to perform the deconvolution. To simulate the effect of 
mis-estimated time-constants, we introduce both a random scatter 
(with RMS = 10 per cent across detectors) and a systematic offset 
(t± identically offset by 10 per cent for all detectors) in the time- 
constants when creating the simulated data. 

• Pointing errors. We simulate the effects of both a random jitter 
and a slow wander in the overall pointing of the telescope by intro- 
ducing a random scatter uncorrelated between time samples (with 
RMS 30 arcsec) and an overall drift in the pointing (1 arcmin drift 
from true pointing over the course of a two-hour observation) when 
creating the simulated time-stream. Once again, the simulated data 
is subsequently analysed assuming perfect pointing registers. 

• Differential transmittance in the HWP. As a simple example 
of a HWP-induced systematic, we have considered a differential 
transmittance of the two linear polarizations by the HWP. Prelimi- 
nary measurements of the QOVER HWPs suggest the level of dif- 
ferential transmittance will be in the region of 1—2 per cent. For 
this work, we consider a 2 per cent differential transmittance in the 
HWP. 

The range of systematic effects we include is not exhaustive. 
In particular, we ignore effects in the HWP, when present, except 
for a mis-estimation of the rotation angle and a differential trans- 
mittance of the two linear polarizations. In addition, we ignore all 
optical effects. In practice, there are many possible HWP-related 
systematic effects which we have not yet considered. In general, 
a thorough analysis of HWP-related systematics require detailed 
physical optics modelling which is beyond the scope of our cur- 
rent analysis. We therefore leave a detailed investigation of HWP- 
related systematics to future work and simply urge the reader to 
bear in mind that where our analysis has included a HWP, we have, 
in most cases, assumed a perfect one. For other optical effects, we 
note that O'Dea, Challinor & Johnson (2007) have already investi- 
gated some relevant effects using analytic and numerical techniques 
and we are currently adapting their flat-sky numerical analysis to 
work with the full-sky simulations described here. Our conclusions 
on the ability of modulation to mitigate systematic effects associ- 
ated with imperfect optics, which will be based on detailed physical 
optics simulations of the QOVER beams (Johnson et al., in prep), 
will be presented in a future paper. Note that there are important 
instrument-specific issues to consider in such a study to do with 
where the modulation is performed in the instrument. In QOVER, 



the modulation will be provided by a HWP between the horns and 
mirrors and this may lead to a difference in the relative rotation of 
the field directions on the sky and the beam shapes as the HWP ro- 
tates compared to a set-up, as in spider (Crill et al. 2008), where the 
HWP is after (thinking in emission) the beam-defining elements. 

For a QOVER-like receiver, consisting of a HWP followed by a 
polarization analyser (e.g. an orthomode transducer) and detectors, 
we include essentially all relevant systematic effects introduced by 
the receiver. To see this, note that the most general Jones matrix 
describing propagation of the two linear (i.e. x and y) polariza- 
tion states through the polarization analyser is (O'Dea, Challinor 
& Johnson 2007) 

where gi, Q2 and a are small real parameters and ei and £2 are 
small and complex-valued. The detector outputs are proportional 
to the power in the x and y-components of the transmitted field (af- 
ter convolution with the detector response function). To first-order 
in small parameters, only gi, g2 and the real parts of ei and £2 en- 
ter the detected power. In this limit, the perturbed Jones matrix is 
therefore equivalent in terms of the detected power to 

/ 1 + gi \ / cosqi sinai \ 
\ 1 + </2 / \ — sin Q2 cos a?2 J 

where the small angles qi £s Kei and 02 ~ — 5fe2 denote the 
perturbations in the polarization-sensitivity directions introduced 
above, and gi and <?2 are the gain errors. Note that instrument po- 
larization (i.e. leakage from T to detected Q or U) is only generated 
in the receiver through mismatches in the gain at first order, but also 
through I ei 1 2 and | e2 1 2 in an exact calculation. The latter effect is 
not present in the simplified description in terms of offsets in the 
polarization-sensitivity directions. Note also that if we difference 
the outputs of the two detectors in the same pixel, in the presence 
of perturbations ai and Q2 to the polarization sensitivity directions 
we find 

di — d,2 = cos(ai — 02) (Q cos[2((f> — a)] 

+Usm[2(<f>-a)]), (13) 

where a = (ai + 02) /2. This is equivalent to a common rotation 
of the pair by a and a decrease in the polarization efficiency to 

cos(ai — Q2). 



3.7 Polarization modulation 

For our reference simulation, and for each of the systematic ef- 
fects listed above, we simulate the experiment using three different 
strategies for modulating the polarization signal. Firstly, we con- 
sider the case where no explicit modulation of the polarization sig- 
nal is performed - in this case, the only modulation achieved is via 
telescope scanning and the relatively small amount of sky rotation 
provided by the current QOVER observing strategies. In addition, 
we also consider the addition of a half-wave plate, either contin- 
uously rotating or "stepped", placed in front of the focal plane. A 
half-wave plate modulates the polarization signal such that the out- 
put of a single detector (in the detector's local polarization coordi- 
nate frame) is 

di = \[T + Q cos(40O + U sin(40i)] , (14) 
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linear polarizations. We model this effect using a non-ideal Jones 
matrix for the HWP of the form, 



k (Hz) 

Figure 5. Frequency-space representation of polarization modulation with 
a continuously rotating HWP. The plotted power spectrum is that for a sin- 
gle azimuth scan from one of our signal-only simulations with the HWP 
continuously rotating at /a =3 Hz. The power in the frequency range 
0.01 < k < 1 Hz is the unmodulated temperature signal from sky multi- 
poles in the range 20 < i < 2000. The power spike at 4f x = 12 Hz is the 
modulated polarized signal. The dashed (dotted) line shows the expected 
temperature (polarization) signal band (with arbitrary normalisation) ap- 
propriate for the scan speed, modulation frequency and beam size we have 
used. The residual power between 1 and 6 Hz and on either side of the 
polarization power spike is due to pixelisation effects. 



where, here, <j>i is the angle between the detector's local polariza- 
tion frame and the principal axes of the wave plate. 

For a continuously rotating HWP, the polarized-sky signal is 
thus modulated at 4/a where fx is the rotation frequency of the 
HWP. As well as allowing all three Stokes parameters to be mea- 
sured from a single detector, modulation with a continuously rotat- 
ing HWP (which we term "fast" modulation in this paper) moves 
the polarization sky-signal to higher frequency and thus away from 
any low-frequency 1// detector noise that may be present; see 
Fig. 5. (Note that the temperature signal is not modulated and one 
needs to rely on telescope scanning and analysis techniques to mit- 
igate 1/f noise in T.) This ability to mitigate the effect of 1// 
noise on the polarization signal is the prime motivation for includ- 
ing a continuously rotating HWP in a CMB polarization exper- 
iment 4 Systematic effects that generate an apparent polarization 
signal that is not modulated at 4/a can also be mitigated almost 
completely with fast modulation. Most notably, instrument polar- 
ization generated in the receiver will not produce a spurious po- 
larization signal in the recovered maps unless the gain and time- 
constant mismatches vary sufficiently rapidly (~ 4/a) to move the 
scan-modulated temperature leakage up into the polarization signal 
band. 

As mentioned in the previous section, as an example of a 
HWP-induced systematic, we have considered the case of a 2 per 
cent differential transmittance by the HWP of the two incoming 



4 The ability to measure all three Stokes parameters from a single detector 
has also been suggested as motivation for including a modulation scheme in 
CMB polarization experiments. However, we will argue later in Section 6.2 
that, at least for ground-based experiments, an analysis based on extracting 
all three Stokes parameters from individual detectors in isolation using a 
real-space demodulation technique may be a poor choice. 



J = 



1 

-(1 + 8) 



(15) 



where 8 describes the level of differential transmittance. Propagat- 
ing through to detected power, for the difference in output of the 
two detectors within a pixel, we find 



tii — d-2 



1 + S + — 

4 



(Qcos(4<?(>) + Usm(4</>)) 



r 521 

5 + Y 



7cos(20) + —Q. 



(16) 



Note that, in this expression, both the HWP angle, <f>, and the Stokes 
parameters are defined in the pixel basis. The first term in equation 
(16) is the ideal detector-differenced signal but mis-calibrated by a 
factor, 8 + 8 2 /4. For reasonable values of 8, this mis-calibration is 
small (< 2 per cent) and, in any case, is easily dealt with during 
a likelihood analysis of the power spectra. The potential problem 
term is the middle term which contains the total intensity signal 
modulated at 2 f\. Note that there will be contributions to this term 
from the CMB monopole, dipole and the atmosphere, which, for 
our simulations, we have taken to be 5 per cent emissive. Even for 
small values of 8 therefore, these HWP-synchronous signals will 
completely dominate the raw detector data and need to be removed 
from the data prior to the map-making step. 

With a HWP operating in "stepped" mode where the angle of 
the HWP is changed at regular time intervals (e.g. at the end of each 
scan), the gains are less clear. The polarization sky signal is not 
shifted to higher frequencies so 1// detector noise can only be dealt 
with by fast scanning. What stepping the waveplate can potentially 
do is to increase the range of polarization sensitivity directions with 
which a given pair of detectors samples any pixel on the sky. This 
has two important effects: (i) it reduces the correlations between 
the errors in the reconstructed Q and U Stokes parameters in each 
sky pixel; and (ii) it can mitigate somewhat those systematic effects 
that do not transform like a true polarization under rotation of the 
waveplate. Of course, one of the strongest motivations for stepped, 
slow modulation is the avoidance of systematic effects associated 
with the continuous rotation of the HWP. If these effects are suffi- 
ciently well understood, then the resulting spurious signals can be 
rejected during analysis. However, if they are not well understood, 
a stepped HWP, while not as effective in mitigating systematics, 
may well be the preferred option. 

Note that for a perfect optical system, rotation of the wave- 
plate is equivalent to rotation (by twice the angle) of the instru- 
ment. However, this need not hold with imperfect optics. For ex- 
ample, suppose the beam patterns for the two polarizations of a 
given feedhorn are purely co-polar (i.e. the polarization sensitivity 
directions are "constant" across the beams and orthogonal), but the 
beam shapes are orthogonal ellipses. This set-up generates instru- 
ment polarization with the result that a temperature distribution that 
is locally quadrupolar on the sky will generate spurious polariza- 
tion that transforms like true sky polarization under rotation of the 
instrument (Hu, Hedman & Zaldariagga 2003; O'Dea, Challinor 
& Johnson 2007). However, for an optical arrangement like that in 
spider, where the HWP is after (in emission) the beam-defining op- 
tics, as the HWP rotates the polarization directions rotate on the sky 
but the beam shapes remain fixed. The spurious polarization from 
the mis-match of beam shapes is then constant as the HWP rotates 
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for any temperature distribution on the sky, and so the quadrupolar 
temperature leakage can be reduced. 

In our analysis, in addition to simulations without explicit 
modulation we have also simulated an experiment with a HWP con- 
tinuously rotating at 3 Hz (thus modulating the polarization signal 
at 12 Hz) and an experiment where a HWP is stepped (by 20°) at 
the end of each azimuth scan (for the scan strategy and scan speed 
we are using, this corresponds to stepping the HWP roughly every 
~90s). 

We end this section with a comment on the ability of a con- 
tinuously rotating HWP to mitigate un-polarized 1// atmospheric 
noise. The polarization signal band is still, of course, moved to 
higher frequency and thus away from the 1 // noise but the atmo- 
spheric 1/ / noise is so strong that extremely rapid HWP rotation 
would be required to move the polarization band far enough into 
the tail of the 1/f spectrum. Such rapid rotation is not an option in 
practice as it would introduce its own systematic effects (e.g. exces- 
sive heat generation). This is the basis of our argument mentioned 
above that extracting all three Stokes parameters from a single de- 
tector may be a poor choice of analysis technique. However, since 
the 1 // atmospheric noise is un-polarized, it can be removed com- 
pletely by combining data from multiple detectors. We revisit this 
issue again with simulations in Section 6.2. Finally, we note that if 
the atmosphere does contain a polarized 1/f component, then we 
expect that this will not be mitigated by modulation - the polarized 
atmosphere would be modulated in the same way as the sky signal 
and would shift up in frequency accordingly. 



4 ANALYSIS OF SIMULATED DATA 

For our reference simulation, and for each of the systematic ef- 
fects and modulation strategies described in the previous section, 
we have created a suite of 50 signal-only, noise-only and signal- 
plus-noise simulated datasets. Our analysis of the signal-only data 
will be used to investigate potential biases caused by the systemat- 
ics while our signal-plus-noise realisations are used together with 
the noise-only simulations to investigate any degradation of the 
sensitivity of the experiment due to the presence of the systematic 
effects. 

Our analysis of each dataset consists of processing the data 
through the stages of deconvolution for the bolometer response 
function, polarization demodulation and map-making, and finally 
estimation of the E- and B-mode power spectra. For any given sin- 
gle realisation these processes are performed separately for each 
of the four observed QOVER fields - that is, we make maps and 
measure power spectra for each field separately. Since our fields 
are widely separated on the sky, we can treat them as independent 
and combine the power spectra measured from each using a simple 
weighted average to produce a single set of E- and B-mode power 
spectra for each realisation of the experiment. Note that even if our 
fields were not widely separated, our procedure would still be unbi- 
ased (but sub-optimal) and correlations between the fields would be 
automatically taken into account in our error analysis since, for any 
given realisation, the input signal for all four fields is taken from 
the same simulated CMB sky. 



4.1 Time-stream processing and map-making 

We first deconvolve the time-stream data for the detector response 
in frequency space using the response function of equation (7) and 




k (Hz) 

Figure 6. Power spectra of a six minute segment of time-stream data before 
(upper panel) and after (lower panel) fitting for and removing the HWP- 
synchronous signal. A 2 per cent differential transmittance in the HWP was 
used to create the simulated data. The resulting spurious signal appears at 
2 fx = 6 Hz in the upper panel and has a peak power of ~ 10 10 in the units 
plotted. No obvious residuals are apparent in the lower panel after applying 
our procedure for removing the contamination. 



using the nominal time-constants in all cases. Once this is done, 
the data from detectors within each pixel are differenced in order 
to remove both the CMB temperature signal and the correlated 1 // 
component of the atmospheric noise. For the case where the 1/f 
atmosphere is completely correlated between the two detectors and 
in the absence of instrumental polarization and/or calibration sys- 
tematics, this process will remove the CMB T signal and the 1/f 
atmosphere completely. 

For the case where we have simulated the effect of a differen- 
tial transmittance in the HWP, a further time-stream processing step 
is required at this point to fit for and remove the HWP-synchronous 
signals from the time-stream. To do this, we have implemented a 
simple iterative least-squares estimator to fit, in turn, for the ampli- 
tudes of both a cosine and sine term at the second harmonic of the 
HWP-rotation frequency, fx. For our simulations containing both 
signal and noise, the accuracy with which we are able to remove 
the HWP-synchronous signals is determined by the noise level in 
the data. A demonstration of the performance of this procedure is 
given in Fig. 6 where we plot the power spectra of six minutes of 
simulated time-stream data (containing both signal and noise) be- 
fore and after the removal of the HWP-synchronous signal. 

After detector differencing, and the removal of the HWP- 
synchronous signals if present, the resulting differenced time 
stream is then a pure polarization signal: 

d i = Qcos(2c/> i ) + Usm(2c/> i ), (17) 

where again, the angle <j>i is, in the most general case, a combina- 
tion of detector orientation, sky crossing angle, telescope boresight 
rotation and the orientation of the HWP if present. In order to con- 
struct maps of the Q and U Stokes parameters, these quantities need 
to be decorrelated from the differenced time-stream using multiple 
observations of the same region of sky taken with different values 
of (f>i. For an experiment which does not continuously modulate the 
polarization signal, the Q and U signals have to be demodulated as 
part of the map-making step. Note however that for an experiment 
where the polarization signal is continuously modulated, there are 
a number of alternative techniques to demodulate Q and U at the 
time-stream level. In separate work, one of us has compared the 
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performance of a number of such demodulation schemes and our 
results will be presented in a forthcoming paper (Brown, in prep.). 
For the purposes of our current analysis however, we have applied 
the same map-based demodulation scheme to all three experiments 
which we have simulated. In this scheme, once the time-streams 
from each detector pair have been differenced, Q and U maps are 
constructed as 



(cos 2 (200) 
(cos(20i) sin(2 < /> i )> 

(cos(2(f>i)di) 
(sin(2(f>i)di) 



(cos(2^i) sin(20i)} 
<sm 2 (2<^)> 



(18) 



where the angled brackets denote an average over all data falling 
within each map pixel. For the work presented here, this averaging 
is performed using the data from all detector pairs in one operation. 
One could alternatively make maps per detector pair which could 
then be co-added later. For the case where the noise properties of 
each detector pair are similar, the two approaches should be equiv- 
alent. Note that the map-maker we use for all of our analyses is 
a naive one - that is, we use simple binning to implement equa- 
tion (18). There are, of course, more optimal techniques available 
(e.g. Sutton et al. 2009 and references therein) which would out- 
perform a naive map-maker in the presence of non-white noise. 
However, for our purposes, where we wish to investigate the im- 
pact of modulation in isolation, it is more appropriate to apply the 
same nai ve map-making technique to all of our simulated data. We 
can then be sure that any improvement we see in results from our 
simulations including modulation are solely due to the modulation 
scheme employed. 



4.2 Power spectrum estimation 

We measure E- and B-mode power spectra from each of our recon- 
structed maps using the "pure" pseudo-CV method of Smith (2006). 
We will not describe the method in detail here and refer the inter- 
ested reader to Smith (2006) and Smith & Zaldarriaga (2007) for 
further details. Here, we simply note that the "pure" pseudo-CV 
framework satisfies most of the requirements of a power spectrum 
estimator for a mega-pixel CMB polarization experiment with com- 
plicated noise properties targeted at constraining B-modes: it is, 
just like normal pseudo-CV, a fast estimator scaling as N*(£ where 
iVpix is the number of map pixels (as opposed to a maximum like- 
lihood estimator which scales as iVp ix ); it is a Monte-Carlo based 
estimator relying on simulations of the noise properties of the ex- 
periment to remove the noise bias and estimate band-power errors 
and covariances - it is thus naturally suited to experiments with 
complicated noise properties for which approximations to the noise 
cannot be made; and it is near-optimal in the sense that it elimi- 
nates excess sample variance from E — > B mixing due to ambigu- 
ous modes which result from incomplete sky observations (Lewis, 
Challinor & Turok 2002; Bunn 2002), and which renders simple 
pseudo-C^ techniques unsuitable for small survey areas (Challinor 
& Chon 2005). 



4.2. 1 Power spectrum weight functions 

With normal pseudo-CV estimators, one multiplies the data with a 
function W(n) that is chosen heuristically and apodizes the edge 
of the survey (to reduce mode-coupling effects). For example, if 
one is signal dominated, uniform weighting (plus apodization) is 



a reasonable choice, whereas an inverse- variance weight is a good 
choice in the noise-dominated regime. Similar reasoning applies for 
the pure pseudo-C^ technique, but here one weights the spherical 
harmonic functions rather than the data themselves. To see this, 
compare the definition of the ordinary and pure pseudo harmonic 
B-modes: 



(1-2)1 f ,2- 



H + (n)W(n)mY e * m (n) 



-U-(n)W(n)dW; m (n) 



(19) 



~B pure 



a 



Cm 



II_(n)83(w(n)Y/ m (n)) 



(20) 



Here, n± (h) = (Q ± ill) (n) is the complex polarization and 9, 9 
are the spin raising and lowering operators defined in Zaldarriaga & 
Seljak (1997). If W(n) is chosen so that it vanishes along with its 
first derivative on the survey boundary, then the afj^ urc couple only 
to B-modes and the excess sample variance due to E-B mixing is 
eliminated. The action of 9, 9 on the spin spherical harmonics is 
simply to convert between different spin-harmonics but their action 
on a general weight function is non-trivial for W(fi) defined on an 
irregular pixelisation such as healpix 5 . To get around this prob- 
lem, we choose to calculate the derivatives of W(fi) in the flat-sky 
approximation where the differential operators reduce to 



9 « -(9a: + id y ), 
9 ~ —(d x — id y ). 



(21) 
(22) 



The derivatives are then trivially calculated on a regular Cartesian 
grid using finite differencing (Smith & Zaldarriaga 2007). 

The most optimal weighting scheme for a pseudo-C"^ anal- 
ysis involves different weight functions for each C'e band-power 
according to the signal-to-noise level expected in that band. How- 
ever, this is a costly solution (requiring 3-/V ban d spherical harmonic 
transforms) and the indications are, from some restricted tests that 
we have carried out, that the improvement in resulting error-bars is 
small, at least for the specific noise properties of our simulations. 
For the analysis presented here, we have therefore chosen a sim- 
pler scheme whereby we have used a uniform weight, appropriately 
apodized at the boundaries for the entire I range for B-modes and 
for I ^ 200 for B-modes. For I > 200 our simulated experiment 
is completely noise dominated for a measurement of B-modes and 
so here we use an inverse- variance weighting, again, appropriately 
apodized at the boundaries. For simplicity, we have approximated 
the boundary of the map as a circle of radius 11° . Note that restrict- 
ing our power spectrum analysis to this central region of our maps 
means we are effectively using only ~ 70 per cent of the available 
data. To calculate the derivatives of the weight functions, we project 
our weight maps (defined in healpix) onto a Cartesian grid using 
a gnomonic projection. Once the derivatives of the weight maps 
have been constructed on the grid using finite differencing, they are 
transformed back to the original healpix grid. An example of the 



5 One possibility that we have yet to explore is performing the derivatives 
directly in spherical-harmonic space. Since Win) is typically smooth, its 
spherical transform will be band-limited and straightforward to handle. 
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Figure 7. Inverse- variance weight functions used for power spectrum estimation for one of the southern fields. For the noise properties of our simulated data, 
the hit-map shown in the top left panel closely approximates the inverse- variance map. This map is heavily smoothed and apodized at the boundary of the map 
to produce the spin-0 weight function shown in the top right panel. The spin-1 and spin-2 weight functions, QW and 99iy, are shown (as vector fields) in the 
bottom left and right-hand panels respectively. 



inverse variance weight maps we have used and the resulting spin-1 
and spin-2 weight functions, dW and 33W\ for one of our fields 
are shown in Fig. 7. 



5 RESULTS FROM SIMULATIONS 

The map-making and power spectrum estimation procedures de- 
scribed above have been applied to each of our simulated datasets 
treating each of our four observing fields independently. For any 
given simulation set, we have 50 Monte-Carlo simulations so we 
can estimate the uncertainties on the band-powers measured from 
each field. For each realisation, we can then combine our measure- 
ments from the four fields using inverse-variance weights to pro- 
duce a final single estimate of the power spectrum for each real- 
isation. When presenting our results below, in all cases, we plot 
the mean of these final estimates. For our reference simulation, and 
for the 1// noise systematics, the error-bars plotted are calculated 
from the scatter among the realisations and are those appropriate 
for a single realisation. When investigating the B-mode bias from 
systematics, we plot the results from signal-only simulations and 
the error-bars plotted are the standard error on the mean. For some 
of the noise-related systematics, we will also examine the impact of 
modulation in the map domain where the effects are already clear. 



5.1 Reference simulation 

To provide a reference for the results which follow, in Figs. 8 and 9 
we show the results from our suite of simulations with no explicit 
modulation and with no systematics included. In Fig. 8 we plot ex- 
amples of the reconstructed noise-only and signal-plus-noise T and 
U Stokes parameter maps (the reconstructed Q maps - not shown 
- are qualitatively similar to U). The raw collecting power of an 
experiment like Clover is apparent from the top two panels in this 
figure. Although the noise T map shown in the top-left panel is 
clearly dominated by striping due to the correlated noise from the 
atmosphere, only signal is apparent in the signal-plus-noise map 
shown in the top-right panel. (In fact, the noise contribution to the 
T signal-plus-noise map is significant, particularly on large scales, 
and so would need to be accounted for when measuring the temper- 
ature power spectrum.) Conversely, the U noise map is dominated 
by white noise; the correlated component of the atmosphere has 
been removed completely from the polarization time-streams (as 
has the T sky signal) by differencing detector pairs before map- 
making. 

Figure 9 shows the mean recovered E- and B-mode power 
spectra from our reference simulations for the case of no explicit 
modulation. Here, we see that our analysis is unbiased and recovers 
the input polarization power spectra correctly. For an input tensor- 
to-scalar ratio of r = 0.026, we recover a detection of B-modes in 
excess of the lensing signal of 1.54er. (We argue in Section 6.3 that 
this is an under-estimate of the detection significance by around 
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Figure 8. Sample maps constructed from simulated time-stream containing noise only (left panels) and both signal and noise (right panels). Temperature maps 
are shown in the top panel and (/-polarization maps are shown in the bottom panels. These maps are for one of our reference simulations with no explicit 
modulation scheme and no systematics included. Note the striping in the noise-only T map which is completely absent from the U maps due to differencing 
of detector pairs before map-making. 



CD 




multipole, I 

Figure 9. Mean recovered i?-mode (top) and B-mode (bottom) power spectra for the reference simulations without explicit modulation. The errors plotted are 
those appropriate for a single realisation. The input CMB power spectra used to create the signal component of the simulations are shown as the red curves. In 
the bottom panel the total B-mode input signal (including lensing) for a tensor-to-scalar ratio of r = 0.026 is shown as the red curve and the .B-mode signal 
due to lensing alone is shown as the dashed curve. The i < 200 multipole range is shown in detail in the inset plots. 
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10 per cent due to our ignoring small anti-correlations between the 
errors in adjacent band-powers.) 

The corresponding plots for the stepped and continuously ro- 
tating HWP are very similar apart from the reconstructed polariza- 
tion maps at the very edges of the fields where a modulation scheme 
increases the ability to decompose into the Q and U Stokes parame- 
ters. Since the edges of the field are excluded in our power spectrum 
analysis in any case (see Fig. 7), we find that the performance (in 
terms of Ce errors) of all three types of experiments which we have 
considered is qualitatively the same in the absence of systematic 
effects. Note that for all the systematics we have considered, the 
effects on the recovery of the B-mode spectrum is negligible and 
so, in the following sections, we plot only the recovered B-mode 
power spectra which are the main focus of this paper. 



5.2 1/f detector noise 

Figure 10 shows the recovered noise-only maps from a simulation 
containing a correlated 1/f detector noise component with a knee 
frequency of /k nee = 0.1 Hz. In this figure, we have plotted the 
noise-only maps from the three types of experiment we have con- 
sidered: no modulation; a HWP which is stepped by 20° at the end 
of each azimuth scan; and a HWP continuously rotating at 3 Hz. 
The impact of modulation on 1/f detector noise is clear from this 
plot - as described in Section 3.7, the continuously rotating HWP 
shifts the polarization band in frequency away from the 1/f detec- 
tor noise leaving only white noise in the resulting map. A stepped 
HWP, on the other hand, does not mitigate 1/f detector noise in 
this way and so noise striping is apparent in the middle panel of 
Fig. 10. 

The B-mode power spectra measured from our signal-plus- 
noise simulations including 1/f detector noise are shown in Fig. 1 1 
(again for /knee = 0.1 Hz) where we show the results from all 
three types of experiment. Examination of the figure suggests that 
the presence of a 1/f component in the detector noise leads to 
a significant degradation in the ability of the unmodulated and 
stepped-HWP experiments to recover the input B-mode signal. 
This degradation happens at all multipoles but is particularly acute 
on the largest scales (£ < 200) where the primordial B-mode sig- 
nal resides. The marginal detection of primordial B-modes (for 
r — 0.026) which we saw in our reference simulation (Fig. 9) is 
now completely destroyed by the presence of the 1/f correlated 
detector noise. Furthermore, note that our analysis of the simulated 
data sets is optimistic in the sense that we have assumed that any 
correlated noise can be modelled accurately - that is, our noise-only 
simulations, which we use to measure the noise bias, are generated 
from the same model noise power spectrum used to generate the 
noise component in our signal-plus-noise simulations. (This is the 
reason that our recovered spectra are unbiased.) However, for the 
analysis of a real experiment, the noise properties need to be mea- 
sured from the real data and there are uncertainties and approxi- 
mations inherent in this process. Any correlated detector noise en- 
countered in a real experiment is unlikely to be understood to the 
level which we have assumed in our analysis and so will likely re- 
sult not only in the increased uncertainties we have demonstrated 
here but also in a biased result at low multipoles. Estimating cross- 
spectra between maps made from subsets of detectors for which the 
1 // detector noise is measured to be uncorrelated is a simple way 
to avoid this noise bias issue, at the expense of a small increase in 
the error-bars (Hinshaw et al. 2003). 



The results from the simulations where we have continuously 
modulated the polarization signal recover the input B-mode sig- 
nal to the same precision that we saw with our reference simula- 
tion - our marginal detection of r = 0.026 is retained even in the 
presence of the correlated detector noise. In Section 6.1 and Ta- 
ble 1 we show quantitatively that, for a detector knee frequency 
of 0.1 Hz, the significance with which the continuously modulated 
experiment detects the primordial B-mode signal is roughly twice 
that found for the un-modulated and stepped-HWP experiments. 
Also detailed in Table 1 are the results from our 1/f noise sim- 
ulations with knee frequencies of 0.05 and 0.01 Hz. We see, as 
expected, that the impact of fast modulation is less for a lower 
knee frequency — for / knoe = 0.05, rapid modulation still signif- 
icantly out-performs the un-modulated and stepped-HWP experi- 
ments while for /knee = 0.01 Hz, there is essentially no difference 
between the performance of the three types of experiment. 

Note that, in the case of rapid modulation, because the polar- 
ization signal is moved completely away from the 1/f frequency 
regime, the recovered spectra should be immune to the issues of 
mis-estimation or poor knowledge of the noise power spectrum at 
low frequencies mentioned above. Although detector 1/f noise 
can be mitigated by other methods (e.g. using a more sophisti- 
cated map-maker; Sutton et al. 2009), these usually require accurate 
knowledge of the low-frequency noise spectrum unlike the hard- 
ware approach of fast modulation. 

5.3 Polarized atmospheric 1/f 

In contrast to the addition of 1// detector noise, which can be suc- 
cessfully dealt with by rapid modulation, all three types of exper- 
iment are degraded similarly by polarized low-frequency noise in 
the atmosphere. In particular, the errors at low multipoles are in- 
flated by a large factor since the large amount of polarized atmo- 
sphere which we have input to the simulations swamps the input B- 
mode signal for r = 0.026. We stress that the levels of polarized 
atmosphere we have used in these simulations were deliberately 
chosen to demonstrate the point that modulation does not help and 
the levels are certainly pessimistic. 

5.4 Calibration errors 

The power spectra recovered from signal-only simulations where 
we introduced random gain errors (constant in time) across the fo- 
cal plane, or 1 per cent systematic A/B gain mis-matches between 
the two detectors within each pixel are shown in Fig. 12. In both 
cases, we see a clear bias in the recovered B-mode signal in the 
absence of fast modulation, but the bias is mitigated entirely by the 
presence of a HWP rotating at 3 Hz. The bias is also mitigated to 
some degree by the stepped HWP but not completely. In our sim- 
ulations, the bias is generally larger for the case of random gain 
errors since the variance (across the focal plane) of the gain mis- 
matches is twice as large in the former case. For our simulations 
where we allowed detector gains to drift over the course of a two- 
hour observation, we found a similarly behaved B-mode bias to 
those shown in Fig. 12 but with a smaller magnitude (due to the 
two-hour drifts averaging down over the eight-hour observation). 

A mis-match between the gains of the two detectors within a 
pixel corresponds to a T — > Q leakage in the detector basis. The 
projection of this instrumental polarization onto the sky will there- 
fore be suppressed if a wide range of sensitivity directions <f>i con- 
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Figure 10. Sample noise-only [/-maps from simulations containing a 1// component correlated to the detector noise 

. For display purposes only, the maps have been smoothed with a Gaussian with a FWHM of 7 arcmin. In the case where explicit modulation is either absent 
(left panel) or slow (stepped HWP; middle panel), the 1// noise results in faint residual stripes in the polarization maps. In the case of a continuously rotating 
HWP, the polarization signal is modulated away from the low frequency 1// resulting in white-noise behaviour in the polarization map (right panel). 
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Figure 11. Recovered B-mode power spectra for the simulations including a correlated 1/f component to the detector noise with /knee = 0.1 Hz for no 
modulation (top), a stepped HWP (middle) and a HWP continuously rotating at 3 Hz (bottom). See Table 1 for the significances with which each experiment 
detects the input signals. 



tribute to each sky pixel, as is the case for fast modulation. Note 
that in the case of a stepped HWP, one should be careful to de- 
sign the stepping strategy in such a way that it does not undo some 
of the effect of sky rotation. During our analysis, we have found 
that the performance of a stepped HWP in mitigating systematic 
effects can depend critically on the direction, magnitude and fre- 
quency of the HWP step applied. In fact, for some set-ups we have 



investigated, a stepped HWP actually worsened the performance in 
comparison to the no modulation case due to interactions between 
the scan strategy and HWP stepping strategy. However, the results 
plotted in Fig. 12 for the stepped HWP case are for a HWP step of 
20° between each azimuth scan which is large and frequent enough 
to ensure that such interactions between the stepping strategy and 
the scan strategy are sub-dominant. 
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Figure 12. Mean recovered B-mode power spectra for the simulations including random gain errors across the focal plane (black points) or systematic 1 per 
cent A/B gain mis-matches between detector pairs (blue points) for no modulation (top), a stepped HWP (middle) and a continuously rotating HWP (bottom). 
These spectra are measured from our suite of signal-only simulations. Our simulations containing both signal and noise exhibit the same biased recovery 
for the no-modulation and stepped-HWP cases. The presence of a fast modulation scheme (bottom panel) mitigates entirely the bias caused by these gain 
mis-matches. The standard errors in these mean recovered spectra are smaller than the plotted symbols. 



5.5 Mis-estimated polarization angles 

The next set of systematics we have considered concern a mis- 
estimation of both the detector orientation angles (i.e. the direction 
of linear polarization to which each detector is sensitive to) and, for 
the case where a HWP is employed, a mis-estimation of the HWP 
orientation. We have performed simulations including both a ran- 
dom scatter (with an RMS of 0.5°) and a systematic offset of 0.5° 
in the simulated detector and HWP angles. Note that for the sys- 
tematic offset in the detector angles, the same offset is applied to 
all detectors. For both the detector angles and the HWP, the offset 
introduced corresponds to a systematic error in the estimation of 
the global polarization coordinate frame of the experiment and the 
effects are therefore degenerate. 

For the simulations which included a random scatter in the an- 
gles (both detector angles and HWP orientation), we found neither 
a bias in the recovered B-mode power spectra, nor a degradation 
in the error-bars from the simulations containing both signal and 
noise. Following the discussion in Section 3.6, common errors in 
the detector angles for the pair of detectors in a single focal-plane 
pixel give rise to a rotation of the polarization sensitivity direc- 
tion of the pixel, while differential errors reduce the polarization 
efficiency. For a typical differential scatter of V2 x 0.5°, the re- 
duction in the polarization efficiency (~ 10~ ) is negligible. For 
a given pixel on the sky, the impact of the polarization rotation is 
suppressed by ^iVsampie, where -/V samp i c is the total number of 



samples contributing to that pixel with independent errors in the 
angles. The combination of a large number of detectors and, in the 
case of random HWP angle errors, their assumed short correlation 
time in our simulations renders the effect of small and random scat- 
ter in the angles negligible. 

The results from simulations which included a systematic er- 
ror in the angles are shown in Fig. 13, where we plot the recovered 
B-mode power spectrum from our signal-only simulations. In con- 
trast to the simulations with random scatter, there is a clear mixing 
between E and B due to the systematic mis-calibration of the po- 
larization coordinate reference system of the instrument. A global 
mis-estimation of the polarization direction by an angle ip in the 
reconstructed maps leads to spurious B-modes with 

Cf = sin 2 {2i>)C? « 4^ 2 Cf . (23) 

Note that in Fig. 13, we show the results from the detector angle 
systematic only for the case of the non-modulated experiment but 
the plot is identical for both the stepped and continuously rotating 
HWP - polarization modulation cannot mitigate a mis-calibration 
of detector angles. The fact that the mixing apparent in Fig. 13 is 
greater for the HWP mis-calibration is simply because rotating the 
waveplate by ij) rotates the polarization direction by 2-0. Although 
the spurious B-mode power is most noticeable at high multipoles, 
where £(£ + l)Cf is largest, it is also present on large scales and, 
as is clear from the plot, would bias a measurement of the B-mode 
spectrum at all multipoles. 
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Figure 13. Mean recovered B-mode power spectra for the signal-only 
simulations including mis-estimated detector polarization sensitivity angles 
(black points) and mis-estimated HWP angles (blue points) where the an- 
gles have been systematically offset by 0.5° in both cases. The standard 
errors in these mean recovered spectra are smaller than the plotted symbols. 



5.8 Differential transmittance in the HWP 

The power spectra reconstructed from our simulations which in- 
cluded a 2 per cent differential transmittance in the HWP exhibited 
no degradation in the accuracy of the recovered B-mode signal — 
even the relatively simple recipe which we have used to remove 
the HWP-synchronous signals from the time-stream (see Section 
4.1) appears sufficient to recover the B-mode signal to the same 
accuracy as was seen in our reference simulations. (We quantify 
this statement in the next section where we estimate the detection 
significances with which the different simulations detect the E and 
B-mode signals). As mentioned in Section 3.7, the recovered po- 
larization signal is mis-calibrated by ~ 2 per cent in amplitude (4 
per cent in power). Compared to the random noise however, this 
mis-calibration is a small effect and is easily dealt with during, e.g. 
a cosmological parameter analysis by marginalising over it. 

Note that no prior information on the level of differential trans- 
mittance was used during our analysis of the data. Our technique 
for removing the HWP-synchronous signals is a blind one in this 
sense and should work equally well for other HWP-systematic ef- 
fects that result in spurious signals at harmonics of the HWP rota- 
tion frequency, fx . 



5.6 Mis-estimated time-constants 

The power spectra measured from simulations which included ran- 
dom and systematic errors in the detector time-constants displayed 
neither a bias, nor a degradation in error-bars. This was to be ex- 
pected for the slow scan speed and extremely fast time-constants 
we have considered in this analysis - the response function of the 
CLOVER detectors is effectively phase-preserving with zero atten- 
uation in the frequency band which contains the sky-signal in our 
simulations. Note that this would not necessarily have been the case 
had we considered a much faster scan speed or more rapid polar- 
ization modulation. 



5.7 Pointing errors 

Our analysis of simulations where we have introduced a jitter in the 
pointing and/or an overall wander in the pointing suggest that these 
systematics have only a very small effect on the recovered B-mode 
power spectra, at least for the levels which we have considered (i.e. 
a 30 arcsec random jitter in the pointing and/or an overall wander 
of the pointing by 1 arcmin over the course of a two-hour observa- 
tion). The only observed effect was a slight suppression of the re- 
covered B-mode signal at high multipoles consistent with a slight 
smearing of the effective beam. We note however that the effect 
we observed was extremely small and was only noticeable in our 
signal-only simulations. For our simulations containing noise, the 
effect was completely swamped by the errors due to random noise. 

In principle, pointing errors can also lead to leakage from E to 
B (Hu, Hedman & Zaldariagga 2003; O'Dea, Challinor & Johnson 
2007). In Appendix A we develop a toy-model for the leakage ex- 
pected from random pointing jitter in the case of a scan/modulation 
strategy that produces a uniform spread of polarization sensitiv- 
ity directions in each sky pixel (such as by fast modulating with a 
HWP). The result is a white-noise spectrum of B-modes but, for 
the simulation parameters adopted here, the effect is very small - 
less than 1 per cent of the B-mode power induced by weak gravi- 
tational lensing on large scales. 



6 DISCUSSION 

6.1 Controlling systematics with polarization modulation 

The main goal of the analysis presented in this paper is to demon- 
strate and quantify with simulations the impact of two types of po- 
larization modulation (slow modulation using a stepped HWP and 
rapid modulation with a continuously rotating HWP) on the science 
return of upcoming CMB B-mode experiments in the presence of 
various systematic effects. Although our list of included systemat- 
ics is not an exhaustive one (in particular, we are still investigat- 
ing the case of imperfect optics), we are nevertheless in a position 
to draw some rather general conclusions regarding the usefulness 
of modulation in mitigating systematics. It is, of course, important 
to bear in mind that we have only considered two examples of a 
HWP-related systematic effect (imperfect HWP angles and differ- 
ential transmittance in the plate). The are many more possible ef- 
fects which will need to be well understood and strictly controlled 
if fast polarization modulation with HWPs is to realise its potential. 

( i ) Systematics mitigated by modulation 

• Correlated 1/f detector noise: As expected by the general 
reasoning of Section 3.7, and further borne out by our results from 
simulations, rapid polarization modulation is extremely powerful 
at mitigating a correlated 1/f component in the detector noise. 
Any such 1/f component is not mitigated by a HWP operating 
in stepped mode. 

• Calibration errors: Our results demonstrate that fast modula- 
tion is also useful for mitigating against possible calibration errors 
since it greatly increases the range of directions over which sky po- 
larization is measured in a given pixel. For example, the clear bias 
introduced in our simulations by random gain drifts or systematic 
mis-calibrations between detectors was mitigated entirely by the 
HWP continuously rotating at 3 Hz. This bias was also partly (but 
not completely) mitigated by stepping the HWP by 20° between 
each of our azimuth scans. For some stepping strategies we have in- 
vestigated, the bias actually increased — a poor choice of stepping 
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strategy can actually be worse than having no modulation because 
of interactions between the sky rotation and the HWP orientations. 

( ii) Systematics not mitigated by modulation 

• Polarized 1 // atmosphere: No amount of modulation (rapid 
or slow) will mitigate a polarized 1// component in the atmo- 
sphere. The results from our simulations containing polarized at- 
mosphere are summarised in Table 1. 

• Pointing errors: For our simulations which included point- 
ing errors, the effect on the recovered B-mode power spectra was 
extremely small and was equivalent to a slight smoothing of the 
effective beam. Although the same amount of smoothing was ob- 
served in all the simulations (and so the effect is not mitigated by 
polarization modulation), the effect is negligible for the sensitivities 
and beam sizes considered here. The further leakage of B-mode 
power into B modes due to pointing errors was, as expected (see 
Appendix A) unobservably small in our simulations. 

• Mis-calibration of polarization angles: A polarization mod- 
ulation scheme does not mitigate a systematic error in the calibra- 
tion of the polarization sensitivity directions. Experiments using a 
HWP will require precise and accurate measurements of the HWP 
angle at any given time to avoid the E — > B mixing apparent in 
Fig. 13. 

Our results are in broad agreement with those of a similar 
study by MacTavish et al. (2008) who based their analysis on 
signal-only simulations of the spider experiment. Both MacTavish 
et al. (2008) and this study find that polarization modulation with 
a continuously rotating HWP is extremely effective in mitigating 
the effects of 1// detector noise but that, in the presence of signifi- 
cant 1// noise, the analysis of an experiment where modulation is 
either absent or slow will require near-optimal map-making tech- 
niques. In addition, both studies find that the effect of small and 
random pointing errors on the science return of upcoming B-mode 
experiments is negligible given the experimental sensitivities. The 
two analyses also find that the effect of random errors (with ~ 0.5° 
RMS) in the detector polarization sensitivity angles is negligible 
but that the global polarization coordinate frame of the experiment 
needs to be measured carefully — MacTavish et al. (2008) quote a 
required accuracy of < 0.25° for spider which is consistent with 
the requirement for an unbiased measurement of the B-mode sig- 
nal (for r = 0.026) at I < 300 with QOVER. Finally, both studies 
suggest that, in the absence of fast modulation, relative gain errors 
will also need to be controlled to the < 1 per cent level (although, 
in this paper, we have demonstrated that such gain errors are almost 
entirely mitigated using a fast modulation scheme; see Fig. 12). 

In Table 1, we quantify the impact of modulation on the 1// 
noise systematics we have considered in this work by considering 
the significances with which we detect the B-mode and B-mode 
signals. For comparison, the detection significances in the presence 
of a 2 per cent differential transmittance in the HWP are also pre- 
sented. To calculate the total significance of the detection we com- 
pute the Fisher error on the amplitude of a fiducial spectrum, 
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where cov~} is the inverse of the band-power covariance matrix 
for the given spectrum. For the total significance of a detection of 
E or B-modes, the fiducial band-powers are simply the binned in- 
put power spectra. In order to estimate the significance of a de- 



tection of primordial B-modes, we subtract the lensing contribu- 
tion from the input B-mode power spectra. Because the primordial 
B-mode power spectrum is directly proportional to the tensor-to- 
scalar ratio, r, the significance with which we detect the B-mode 
signal in excess of the lensing signal translates directly to a sig- 
nificance for the detection of our input tensor-to-scalar ratio of 
r = 0.026. When analysing the results of the simulations, we ap- 
proximate cov~J ps S bb i jal since we are unable to estimate the 
off-diagonal elements from our small number of realisations (50) 
in each simulation set. We know from a Fisher-based analysis (see 
Section 6.3), the results of which are also reported in Table 1, that 
neighbouring band-powers on the largest scales are, in fact, ~ 10 
per cent anti-correlated, and the diagonal approximation therefore 
underestimates the detection significance by ~ 10 per cent. For 
consistency, the numbers quoted for the Fisher analysis ignore the 
off-diagonal elements of the covariance matrix. Including the cor- 
relations increases the B-mode significance to 144.6 (from 128.6), 
the total B-mode significance to 11.57 (from 10.24) and the pri- 
mordial B-mode significance to 2.04 (from 1.90). 

In comparing the entries in Table 1 one should keep in mind 
that the significances reported for the simulations are subject to a 
Monte-Carlo error due to the finite number (iV s i m = 50) of sim- 
ulations used to estimate the band-power errors. Approximating 
the band-powers as uncorrelated and Gaussian distributed, the sam- 
pling error in our estimates of the S/N is 

_i j_ r v (ft 4 ! 

S/N V2N~ V V °b 



A(S/N) : 



(25) 



For the reference simulation, this gives an error of 0.15 in the sig- 
nificance of a detection of r and 0.22 in the significance of the total 
B-mode spectrum. The size of these errors likely explain the ap- 
parent anomalies that rotating the HWP degrades the detection of 
r in the reference simulation, and that adding 1// detector noise 
improves the detection of r over the reference simulation for the 
case of a rotating HWP. 

Also included in Table 1 is the performance of our experiment 
as estimated from a set of simple map-based simulations where we 
have injected uniform and isotropic white noise into signal-only T, 
Q and U maps directly. In these simulations, and also for the Fisher 
analysis, the white-noise levels were chosen to match the noise lev- 
els in our main analysis and so they have identical raw sensitiv- 
ity to the time-stream simulations but with perfectly behaved noise 
properties. The broad agreement between our full time-stream sim- 
ulations and these simple map-based simulations suggests that the 
anisotropic noise distribution introduced by the QOVER scan strat- 
egy does not have a large impact on the performance of the experi- 
ment. This agreement also suggests that the slightly poorer perfor- 
mance of the simulations in recovering the r = 0.026 primordial 
B-mode signal as compared to the Fisher predictions is due to the 
sub-optimal performance on large scales of the (pure) pseudo-Ci 
estimator we have used. 



6.2 Importance of combining data from multiple detectors 

For all of our analyses up to this point, in order to remove the cor- 
related 1// atmospheric noise from the polarization analysis, we 
have differenced detector pairs before map-making. However, as 
mentioned in Section 3.7, for the case of a continuously modulated 
experiment, it is possible to measure all three Stokes parameters 
from a single detector in isolation. Here, we argue that this may be 
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Table 1. Detection significances (in units of cr) for our reference simulations, for our simulations with 1// noise systematics and for our simulations with a 
2 per cent differential transmittance in the HWP. Also included for comparison are the predicted detection significances from a Fisher matrix analysis of the 
power spectrum errors (see Section 6.3) and from the simulations containing isotropic and uniform Gaussian noise (see text). The rightmost column displays 
the significance of the detection of the B-mode signal in excess of the lensing signal which corresponds directly to the significance with which each simulation 
detects the input tensor-to-scalar ratio of r = 0.026. 
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a poor choice of analysis technique in the presence of a highly cor- 
related and common-mode systematic such as atmospheric 1/f, at 
least when one employs real-space demodulation techniques such 
as those that we have used in this analysis. The key point to ap- 
preciate here is that, even with a rapid modulation scheme, and 
for an ideal experiment, it is impossible to separate completely 
the temperature and polarization signals in real space using only a 
single detector. In contrast, the technique of detector differencing 
achieves this complete separation of the temperature and polariza- 
tion signals (again for the case of an ideal experiment). Note that 
this is true even in the signal-only case. Consider again the modu- 
lated signal, in the absence of noise, from a single detector sensitive 
to a single polarization: 

d % = [T(<9) + Q{6) cos(2&) + U{6) sin(2<fc)] /2. (26) 

If for each observed data point, di, the true sky signals, T, Q and 
U are different (as is the case for a scanning experiment), there is 
clearly no way to recover the true values of T, Q and U at each 
point in time. The approximation that one must make in order to 
demodulate the data in real space goes to the very heart of map- 
making - that the true continuously varying sky signal can be ap- 
proximated as a pixelised distribution where T, Q and U are taken 
to be constant within each map-pixel. Armed with this assumption, 



6 We note that this is not necessarily true for the case of genuinely band- 
limited temperature and polarization signals when a classical lock-in tech- 
nique (such as that used in the analysis of the MAXIPOL data; Johnson et 
al. 2007) is used to perform the demodulation. We are currently working to 
integrate such a technique into our analysis. 



all three Stokes parameters can be reconstructed from a single de- 
tector time-stream using a generalisation of equation (18): 

/ T \ ( (d t ) \ 

Q = 2M" 1 • (cos(20i)di> . (27) 
V U J V <sin(2&)di) / 

where the decorrelation matrix, M is now given by 

/ 1 (cos(20i)> (sin(20 i )> \ 

M = I (cos(20»)) (cos 2 ^)) (cos(2&) sin(20i)> .(28) 
\(sin(20i)) (cos(20i) sin(2</>i)) (sin 2 (2<^)} / 

If, on the other hand, the data from different detectors are combined 
(e.g. when detector differencing is used), the situation is different - 
because the two detectors within a pixel observe exactly the same 
un-polarized component of the sky signal at exactly the same time, 
differencing the detectors removes the T signal completely without 
any assumptions regarding the scale over which the true sky signal 
is constant. In this case, the decorrelation of Q and U using equa- 
tion (18) still requires an assumption regarding the constancy of the 
Q and U signals over the scale of a map-pixel but now the much 
larger temperature signal has been removed from the polarization 
analysis completely. 

Now, if in addition to the sky signal, we have a common-mode 
time- varying systematic such as an un-polarized 1/f component in 
the atmosphere, this contaminant will again be removed entirely 
with the detector differencing technique (as long as it is completely 
correlated between the two detectors) whilst it will introduce a fur- 
ther approximation into any attempt to decorrelate all three Stokes 
parameters from a single detector time-stream using equation (27). 

To illustrate this point, and to stress the importance of combin- 
ing data from multiple detectors, we have re-analysed the simulated 
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Figure 14. Recovered noise only U -polarization maps from one of our ref- 
erence simulations with continuous modulation. The map on the left is re- 
constructed from demodulated single detector "pure" U time streams and 
has not used information from multiple detectors to separate the T and po- 
larization signals. The map on the right is made from demodulated detector- 
pair "pure" U time-streams and explicitly combines information from the 
two detectors within each pixel to separate T from Q/U. Although striping 
is absent from both maps, the white-noise level in the detector-differenced 
map is reduced compared to that made using the non-differencing analysis. 



data from our set of continuously modulated reference simulations 
but now we perform the demodulation at the time-stream level for 
either single detectors or single detector-pairs in isolation. Firstly, 
we demodulate each detector time-stream individually using equa- 
tion (27) but now the averaging is performed over short segments in 
time rather than over all data falling within a map-pixel. This pro- 
cedure results in "pure" T, Q and U time-streams for each detector 
but at a reduced data rate determined by the number of data sam- 
ples over which the averaging of equation (27) is performed. For 
our second analysis, we first difference each detector pair and then 
demodulate the differenced time-streams using equation (18), again 
applied over short segments of time, resulting in "pure" Q and U 
time-streams for each detector pair, once again at a reduced data 
rate. Maps of the Q and U Stokes parameters are then constructed 
by simple binning of the demodulated Q and U time-streams from 
all detectors or detector pairs. In the case where detector pairs are 
differenced, we are explicitly combining information from multi- 
ple detectors to separate the temperature and polarization signals 
whilst when we do not difference, we are attempting to separate 
the T and Q/U signals present in detector time-stream in isolation. 

The results of these tests are shown in Figs. 14 and 15. Fig- 
ure 14 shows the noise-only [/-polarization maps recovered using 
the two different analysis techniques. Although striping from the 
atmospheric fluctuations is not present in either of the maps, the 
extra uncertainty introduced when one attempts to separate the T 
and Q/U signals from individual detectors in isolation clearly re- 
sults in an increased white noise level in the polarization maps. This 
results in a degradation factor of ~ 2 in the resulting measurements 
of the B-mode power spectrum on all angular scales (Fig. 15) with 
a corresponding degradation in the detection significances for both 
a measurement of the total B-mode signal and for a detection of 
r — 0.026 (Table 2). These results are in excellent agreement with 
those of Sutton et al. (2009) who found it necessary to apply opti- 
mal mapping techniques to rapidly modulated single-detector time- 
streams in the presence of 1// atmospheric noise. 

We emphasise that the results presented in Fig. 15 and in Ta- 
ble 2 for the case where we have analysed single detectors in isola- 
tion would likely be improved if the Fourier domain filtering used 




1000 



multipole, t 



Figure 15. Comparison between the B-mode power spectra recovered us- 
ing an analysis based on differencing detector pairs (black points) and one 
based on demodulating each detector individually (blue points). In the pres- 
ence of a time-varying common-mode systematic, such as the 1/f atmo- 
spheric noise we have considered here, the analysis based on detector dif- 
ferencing is far superior to the analysis based on demodulating each detector 
individually. 



Table 2. Detection significances (in units of a) from the analysis of identical 
simulated data with a HWP rotating continuously at 3 Hz. The first analy- 
sis is based on detector differencing, the second based on demodulation of 
individual detectors in isolation. 



Analysis 



i?-mode B-mode 



: 0.026 



Detector differencing 132.9 10.14 1.43 
Demodulation 123.4 5.11 0.68 



in lohnson et al. (2007) was implemented. We are currently work- 
ing to integrate this step into our algorithm, and we expect to report 
the subsequent improvement in future publications. 

6.3 Comparison of simulated and predicted qover 
performance 

It is common practice to make predictions of the performance of 
upcoming experiments using a Fisher-matrix analysis which at- 
tempts to predict the achievable errors on, for example, power spec- 
tra or cosmological parameters under some simplifying assump- 
tions. Generally these assumptions will include uniform coverage 
of the observing fields and uncorrelated Gaussian noise resulting 
in an isotropic and uniform white-noise distribution across the ob- 
serving field. In contrast, the work described in this paper has made 
use of a detailed simulation pipeline which we have created for the 
QOVER experiment. Our simulation pipeline includes the QOVER 
focal-plane designs as well as a realistic scan strategy appropri- 
ate for observing the four chosen QOVER fields from the telescope 
site in Chile. In addition we have employed a detailed model of 
the TES detector noise properties and responsivity, and 1/f atmo- 
spheric noise correlated across the focal-plane array. Moreover, our 
errors are calculated using a Monte-Carlo analysis and so should 
automatically include any effects due to correlations between map 
pixels etc. An interesting exercise therefore is to compare the ex- 
pected errors from a Fisher-matrix analysis to those obtained from 
our simulation analysis. 
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The polarization band-power Fisher matrix is (e.g. Tegmark & 
de Oliveira-Costa 2001) 



1 / J 8C 8C 



(29) 



where P and P' are B or B, and b labels the bandpower. Here, C 
is the covariance matrix of the noisy Stokes maps and Cjf are band- 
powers of 1(1 + l)Cf /(2ty). We analyse a single circular field 
with area equal to that retained in the pseudo-CV analysis described 
in Section 4.2, and multiply the Fisher matrix by four to account 
for the number of fields observed (which are thus assumed to be 
fully independent). We ignore inhomogeneity of the noise in the 
maps so that our problem has azimuthal symmetry about the field 
centre. This allows us to work in a basis where the data is Fourier 
transformed in azimuth and the covariance matrix becomes block 
diagonal, thus speeding up the computation of the Fisher matrix 
considerably. The Fisher matrix takes full account of band-power 
correlations (both between b and polarization type) and the effect 
of ambiguities in isolating E and B-modes given the survey geom- 
etry. We deal with power on scales larger than the survey by includ- 
ing a junk band-power for each of E and B whose contribution to 
cov(6p)(6p)/ we remove before computing detection significances. 

The comparison between the predicted and simulated perfor- 
mance of QOVER is shown in Fig. 16. In this plot, we also include 
the predicted B-mode errors from a naive mode-counting argument 
based on the fraction of sky observed, /sky For these estimates, we 
assume independent measurements of the power spectrum in bands 
of width A£ given by 

2 



(AC*) = 



-ACi + NtY 



(30) 
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where Ce is the band-averaged input signal and Ni is the band- 
averaged noise. For uncorrelated and isotropic Gaussian random 



noise, the latter is given by Ne 



w 1 B„ 2 where Be 



exp(—£(£ + 1)ct b /2) is the transform of the beam with o\b 
0B/V81n2 for a beam with FWHM of Ob. The weight w' 1 
fipixfpix, where the pixel noise in the Q and U maps is 
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Here O 2 is the total observed area, t ^ B is the total observation time 
and Opi x is the pixel size. In equation (31), we have used NET/\/2 
to account for the fact that a single measurement of Q or U requires 
a measurement from two detectors (or, alternatively, two measure- 
ments from a single detector) and we use iVdet /4 as the effective 
number of Q and U detectors. 

Over most of the i range, the agreement between the Fisher 
matrix predictions and the simulated performance is rather good 
- the only significant discrepancy is for the lowest band-power 
where the simulations fail to match the predicted Fisher error. This 
is almost certainly due to the relatively poor performance of our 
power spectrum estimator on the very largest scales where pseudo- 
Ci techniques are known to be sub-optimal (compared to, for ex- 
ample, a maximum likelihood analysis). In terms of a detection of 
the total B-mode signal, the Fisher analysis predicts a detection for 
QOVER of ~ 12. Oct. For comparison, the naive / s k y analysis pre- 
dicts a 12.4ct detection. For our assumed tensor-to-scalar ratio, the 
Fisher matrix analysis predicts r — 0.026 ± 0.013 (a 2.04ct de- 
tection) and the naive analysis yields r — 0.026 ± 0.011 (a 2.39ct 
detection). Comparing to the detection significances quoted in Ta- 
ble 1, we see that the detections recovered from the simulations fail 
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Figure 16. Comparison between the predicted performance of QOVER as 
calculated using a Fisher-matrix analysis and the simulated performance 
from our Monte-Carlo pipeline (for our reference simulation). Also shown 
for comparison are the errors predicted from a naive f s k y analysis. 



to match these numbers. For the case of the total B-mode ampli- 
tude, this discrepancy is entirely due to the fact that we are unable 
to measure and include in our analysis (anti-)correlations between 
the band-powers measured from our small number (50) of simula- 
tions - when we neglect the correlations in the Fisher matrix anal- 
ysis, the Fisher prediction drops to a 10.24ct detection, in excellent 
agreement with our measured value from simulations. For the pri- 
mordial B-mode signal only, the discrepancy found is also partly 
due to the same effect (ignoring correlations in the Fisher analysis 
reduces the Fisher prediction for primordial B-modes to 1.9ct). As 
mentioned above, we suspect that the additional decrease in sensi- 
tivity to primordial B-modes seen in the simulations is due to the 
slightly sub-optimal performance of our implementation of the pure 
pseudo-C^ estimator on the largest scales. 

We should point out that in this work, we have made no at- 
tempt to optimise the survey strategy in light of recent instrument 
developments. In particular, the survey size we have adopted for 
these simulations was optimised for a measurement of r = 0.01 
with QOVER when the experiment was expected to have twice the 
number of detectors now planned. For the instrument parameters 
we have adopted in this analysis (which are a fair representation 
of the currently envisaged experiment), the optimal survey area for 
a measurement of r = 0.026 would be significantly smaller than 
the ~ 1500 deg 2 we have used here due to the increased noise lev- 
els from the reduced number of detectors. Alternatively, if we had 
assumed a larger input value of r, the optimal survey size would 
increase. Optimisation of both the survey area and the scan strat- 
egy in light of these changes in the instrument design is the sub- 
ject of on-going work. There are, of course, many other sources 
of uncertainty which we have not yet accounted for in our simula- 
tion pipeline and so both the predicted and simulated performance 
numbers should be taken only as guidelines at this time. However, 
it is encouraging that the extra sources of uncertainty which are in- 
cluded in our simulation pipeline (realistic instrument parameters, 
a realistic scan strategy, correlated noise), in addition to any uncer- 
tainties introduced as part of our subsequent analysis of the simu- 
lated data, do not degrade the expected performance of QOVER by 
a large amount. 
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7 CONCLUSIONS 

We have performed a detailed investigation of the ability of both 
slow and fast polarization modulation schemes to mitigate possible 
systematic effects in upcoming CMB polarization experiments, tar- 
geted at measuring the B-mode signature of gravitational waves in 
the early universe. To do this we have used a simulation pipeline 
developed in the context of the Clover experiment, which in- 
cludes realistic instrument and observation parameters as well as 
1// detector noise and 1// atmospheric noise correlated across 
the QOVER focal-plane array. Using this simulation tool, we have 
performed simulations of QOVER operating with no explicit modu- 
lation, with a stepped HWP and with a HWP rotating continuously 
at 3 Hz. We have analysed the resulting time-stream simulations 
using the technique of detector differencing coupled with a naive 
map-making scheme, and finally have reconstructed the E and B- 
mode power spectra using an implementation of the near-optimal 
"pure" pseudo-C^ power spectrum estimator. 

As expected, we find that fast modulation via a continuously 
rotating HWP is extremely powerful in mitigating a correlated 1 // 
component in the detector noise but that a stepped HWP is not. In 
addition, we have demonstrated that a polarized 1// component in 
the atmosphere is not mitigated by any amount of modulation and 
if present, would need to be mitigated in the analysis using a so- 
phisticated map-making technique. We have further verified with 
simulations that fast modulation is very effective in mitigating in- 
strumental polarization that is fixed relative to the instrument ba- 
sis, for example the T — > Q leakage caused by systematic gain 
errors and mis-matches between detectors, in agreement with the 
conclusions of O'Dea, Challinor & Johnson (2007). We have also 
demonstrated that modulation does not mitigate a systematic mis- 
calibration of polarization angles and that these angles will need 
to be measured accurately in order to avoid a systematic leakage 
between E and B-modes. The other systematics which we have 
investigated (pointing errors, mis-estimated time-constants) have a 
negligible impact on the recovered power spectra for the parame- 
ters adopted in our simulations. 

In addition to our investigation of systematic effects, we have 
stressed the importance of combining data from multiple detectors 
and have demonstrated the superior performance of a differencing 
technique as opposed to one based on measuring all three Stokes 
parameters from single detectors in isolation. We suggest that the 
latter technique, although possible in the presence of rapid modu- 
lation, is likely a poor choice of analysis technique, at least in the 
presence of a common-mode systematic effect such as atmospheric 
1// noise. 

Finally, we have compared the simulated performance of the 
QOVER experiment with the expected performance from a simpli- 
fied Fisher-matrix analysis. For all but the very lowest multipoles, 
where the simulations fail to match the Fisher predictions, we find 
excellent agreement between the predicted and simulated perfor- 
mance. In particular, despite the highly anisotropic noise distribu- 
tion present in our simulated maps, our measurement of the total 
B-mode signal matches closely with the Fisher matrix prediction 
(the latter assuming isotropic noise). On the other hand, the mea- 
surement of the large scale B-mode signal (and thus of the tensor- 
to-scalar ratio, r) from the simulations is around 20 per cent worse 
than the Fisher prediction. This is almost certainly due to the sub- 
optimality of our power spectrum analysis on large scales. It is pos- 
sible that the Fisher matrix predictions could be recovered from 
the simulations by using a more optimal weighting scheme in the 



pure pseudo-C^ analysis, or, more likely, by using a maximum- 
likelihood Ce estimator for the low multipoles. 

One important class of systematic effects which we have not 
considered in this paper are those associated with imperfect op- 
tics. Additionally, we have considered only two effects associated 
with an imperfect HWP. The efficacy of fast modulation to mitigate 
systematic effects from imperfect optics, for example instrumental 
polarization due to beam mis-match, is expected to depend crit- 
ically on the optical design (such as the HWP location). We are 
currently working to include such optical systematics in our sim- 
ulation pipeline, along with a more detailed physical model of the 
atmosphere and models of the expected polarized foreground emis- 
sion. In future work, in addition to investigating further systematic 
effects, we will extend our simulations to multi-frequency obser- 
vations and will use these to test alternative foreground removal 
techniques. We will also apply the "destriping" map-making tech- 
nique of Sutton et al. (2009) to our simulations to assess the relative 
merits of destriping in analysis as opposed to a hardware based ap- 
proach for mitigating 1// noise. 
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APPENDIX A: E-B LEAKAGE FROM POINTING JITTER 

For the case of random pointing jitter and a scan/modulation strat- 
egy that ensures every sky pixel is visited with a wide range of po- 
larization sensitivity directions, it is straightforward to estimate the 
spurious B-mode power that is induced from B-modes for a single 
focal-plane pixel. We work in the flat-sky approximation where di- 
rections on the sky are denoted by positions x in the tangent plane. 
A given sky pixel p receives N p hits with polarization sensitivity 
directions {<pi}, where i — 1, . . . , N p , and assume these angles 
are uniformly distributed. If di is the ith differential measurement 
from the detector pair contributing to pixel p, then we can approx- 
imate the recovered Stokes parameters in that pixel (equation 18) 
as 

Q(xp) ~ ^-^d i cos{2(j) i ) 

P i 

U(x p ) w ^- £)di sin(2&). (Al) 



If the ith observation has a pointing error a;, then the signal con- 
tribution to Q(x p ), for example, is 

Q( x p) ~ jj-^2[Q{x P + a l )cos 2 (2(t> i ) 

p i 

+ U(x p + on) sin(20i) cos(2(^)] -(A2) 

Note that here, Q and U are the beam-smoothed fields. Approxi- 
mating the displaced Stokes parameters by a gradient approxima- 
tion, we find 

Q(x p ) « Q{x p ) + a efI • VQ\ Xp + /3 cff ■ VU\ Xp , (A3) 
where the effective angles 



C*cff = 




^eff = j^-^2cos(2(t>i)sm{2ct>i)ai. (A4) 

P i 

Repeating for the recovered U Stokes parameter, we find 

U(x p ) « U(x p ) + 7off ■ VU\ Xp + (3 cS ■ VQU P (A5) 
where 

Teff = Jf J2 sin2 (200"; • (A6) 

P i 

Note that for pointing errors that are not constant in time, Q and 
U are generally displaced differently and there is non-local Q-U 
coupling through /3 cff . We recover the simple map-based model 
of pointing errors used in Hu, Hedman & Zaldariagga (2003) 
and O'Dea, Challinor & Johnson (2007), 

(Q ± ill){x p ) « (Q ± iU)(x p ) + a ■ V(Q ± iU)\ Xp , (A7) 

for the case of constant pointing errors, i.e. a, = a. 

In the limit that the pointing jitter is independent between time 
samples, the pointing errors in the map, a e s, /3 eff and 7 eff are 
independent between pixels. (Note that this will not hold if maps 
are made from more than one focal-plane pixel.) These pointing 
errors are mean zero, and have independent x- and j/-components. 



The non- vanishing correlators are 

(a cff -a eff ) « 3a t 2 od /(2Ag (A8) 

(0 eS -0e«) « a?od/(2iV p ) (A9) 

<7eff7eff> « 3a t 2 od /(2iV p ) (A10) 

<a cff - 7cff ) « a t 2 od /(2iV p ). (All) 



Here, a 2 od = (a, ■ a,) is the variance of the pointing jitter, and 
we have approximated discrete sums over the angles <j>i by their in- 
tegrals. We shall further assume that the number of hits per pixel 
is uniform over the map in which case the map-based pointing er- 
rors have homogeneous statistics through the map. The quantity 
Q tod^pix/A'p will arise in our final result for the spurious B-mode 
power, where f2 p i x is the pixel area. This product is independent 
of pixel area since N p oc f2 p i x . For 30-arcsec jitter that is random 
between 100 Hz samples, a to d \/^a7/N p ~ 3 x 10~ 8 rad 2 for 
an eight-hour observation of a 380-deg 2 field. 

The calculation of the B-mode power due to leakage of E- 
modes follows from equations (A3) and (A5) with the standard 
techniques used, for example, in the calculation of gravitational 
lensing on the CMB (see, e.g. Lewis & Challinor 2006). We shall 
not give the details here but simply note that the final result is a 
white-noise contribution with 
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c b aLA>ix r £de C E -t>„i (A12) 

2N P J 2ir 

where as is the beam size. The integral here approximates the 
mean-squared gradient of the beam-smoothed polarization which, 
for a 5.5-arcmin beam « 1.8 x 10 7 /iK 2 . For our simulation pa- 
rameters we therefore expect Cf < 10~ 8 fiK 2 from pointing jitter 
which is < 1 per cent of the large-scale power induced by gravita- 
tional lensing. 



